Integrated control system for stability control of yaw, roll and lateral motion of a driving vehicle using an integrated sensing system to determine lateral velocity

ABSTRACT

A method of controlling a vehicle includes determining a lateral tire force, a front lateral tire force, a rear lateral tire force, and determining a linear sideslip angle from the front lateral tire force and the rear lateral tire force. The method further includes determining a linear lateral velocity in response to the linear sideslip angle and controlling the vehicle in response to the linear sideslip angle.

TECHNICAL FIELD

The present invention relates generally to a control apparatus for controlling a system of an automotive vehicle in response to sensed dynamic behavior, and more specifically, to a method and apparatus for dynamically controlling a yaw, roll and lateral stability of a vehicle in response to the motion states of the vehicle measured and sensed through an integrated sensing system.

BACKGROUND

Various automotive vehicles have recently begun including vehicle system controls. Vehicle system controls use the sensors to sense a vehicle's dynamics, the driver's intention or even the environmental information around the vehicle; use the electronic control units (ECUs) to process the sensed information; and use the available actuators to conduct suitable control actions requested from the ECUs so as to achieve controlled performance for a driving vehicle. The ongoing goal of the vehicle system controls is to achieve an improved system performance or new functions. Ultimately, the various vehicle control systems need to be coordinately controlled so as to achieve enhanced system level control performances for ride, handling, safety and fuel economy of the vehicle.

One of the vehicle system controls are the various types of vehicle dynamics controls, which use electronically controlled chassis, powertrain or drivetrain subsystems to augment the driver's controllability of the vehicle for ride, handling and safety purposes.

Several existing stability control systems are available. For example, the yaw stability control system (YSC) or a roll stability control system (RSC) is currently equipped on millions of vehicles. The separated stability control functions are the focuses of the existing systems.

It may be desirable to integrate the functions of the various dynamic control systems. One reason that such an integration is desirable is that there is overlap among the operation ranges of the individual control systems. For example, a rollover control system tries to force the vehicle under-steer more during an un-tripped rollover event, while an ESC system might take the vehicle's response under RSC control as an indication of the vehicle's true under-steer to conduct the vehicle's under-steer correction. If there is no integration between RSC and YSC, the control action of those may cancel each out. Hence a successful implementation of RSC and YSC in the same vehicle will need control coordination. One of the enablers of the coordination relies in the accurate discrimination of the features of different unstable vehicle dynamics, which require different control functions.

With current advances in mechatronics, the aforementioned technology enabler is possible. For example, the advanced sensors might be used. Such advanced sensors together with other mechatronics not only help function integration, they also helps improve the individual vehicle stability control performances and provides opportunities for achieving control performances which were previously reserved for spacecraft and aircraft. For example, gyro sensors, previously only used in aircraft and spacecraft, have now been incorporated in various vehicle dynamics controls and the anti-lock brake systems once invented for airplanes are now standard automotive commodities.

With the application of the advanced sensors and mechatronics, superior information about the vehicle's dynamics states can be obtained which can be used to calculate effective feedback and feedforward controls for various functions and their integrations. It is especially effective to use the superior information from the advanced sensors and their innovative sensing algorithms to identify the scenarios, which are otherwise undetectable during a driving involved with complicated road conditions and aggressive driving inputs from the driver.

A typical vehicle stability control and an integrated stability control system sense and control a vehicle's dynamics conditions described in the 3-dimensional space. Those control systems might require the measurements of all or part of the three-dimensional motions which include the rotational motions along the vehicle's roll, pitch, and yaw directions and the translation motion along the vehicle's longitudinal, lateral and vertical directions.

More specifically, in a typical vehicle stability control system such as a yaw stability control system or a roll stability control system, the primary control task is to stabilize the vehicle in yaw or roll directions which will likely involve motions along its roll and yaw directions, and longitudinal and lateral directions. The coupling between different motion directions may not be as strong as in an aircraft or a spacecraft; however, they cannot be neglected in most of the control regions where unstable vehicle dynamics is seen. Unstable vehicle dynamics might involve rolling over or yawing out of the traveling course during aggressive driving when the driver's inputs are well beyond the values allowed by the adhesion limit of the road surface. For example, the excessive steering of a vehicle will likely lead to a unstable yaw and lateral motions, which further cause large rolling motion towards the outside of the turn. If the driver brakes the vehicle during such an excessive steering, the vehicle will also have roll and pitch motions, large load transfer and large lateral discursion. It is desirable for a high performance stability control system to integrate YSC, RSC and lateral stability control (LSC) such that the coupled unstable vehicle dynamics can be stabilized.

Theoretically, if there are sensors which can directly measure the vehicle's roll states, yaw states and lateral states during complicated unstable dynamics, the successful integration would involve: (i) identifying the dominated direction where a primary unstable dynamics can occur; (ii) prioritizing control functions when multiple safety-critical unstable vehicle dynamics can occur; (iii) maximizing control functions when multiple unstable vehicle dynamics are equally critical; (iv) determining transition control actions from one stability control function to another or to a normal operation function. The decision rule for any of the above actions needs close discrimination. Unfortunately, there are no direct measurements of the aforementioned vehicle states. Even the advanced sensors can only provide indirect measurements of the involved vehicle states. Therefore, intelligent sensing algorithms are required.

Therefore, it would be desirable to integrate the YSC, RSC and LSC functions to provide accurate determination of the involved roll, yaw and lateral motions of a vehicle. The controllable variables associated with those motions might be the vehicle's global and relative angular motion variables such as attitudes and the directional motion variables such as the longitudinal and lateral velocities. For instance, in RSC, the relative roll angle between the vehicle body and the road surface is controlled. In LSC, the relative yaw angle between the vehicle's travel direction and the path is controlled. The coupled dynamics among roll, yaw and pitch motions and the elevated (banked or inclined) road caused vehicle motions all need to be differentiated. The driver's intention-based vehicle dynamic behavior will also need to be determined. Due to the fact that the sensors measure the total motion of the vehicle, the sensor measurements contain the total information that includes the values of the driver-induced dynamics, the road geometry-induced dynamics and the dynamics from the gravity together with sensor uncertainties. It is not hard to find that separating the driver-induced dynamics and the road-induced dynamics from the total sensor measurements need intelligent algorithms. Only if such information is determined from the total sensor measurements, the control system can provide adequate control actions in needed situations (driver-induced excessively dynamic or unstable) and generate less occurrences of false control action in unneeded situations such as in certain road geometry-induced dynamics.

In order to achieve the aforementioned separation of the driver maneuver-induced vehicle dynamics information and the road geometry-induced dynamics information from the sensor uncertainties, gravity-induced dynamics terms, a new vehicle sensing technology which contains an inertial measurement unit (IMU) has been pursued at Ford Motor Company. This sensing system is called an Integrated Sensing System (short to ISS) in this invention.

The IMU has been used in inertial navigation system (INS) for aircraft, spacecraft and satellite for decades. Typically an INS system determines the attitude and directional velocity of a flight vehicle through the sensor signals from the IMU sensors and GPS signals, The IMU sensor set includes 3 gyros and 3 linear accelerometers. The INS contains an IMU and a processor unit to compute the navigation solutions necessary for navigation, attitude reference and various other data communication sources. As the same token, the ISS will also be used (but not limited) to determine the vehicle's attitude and directional velocities with the exception that GPS signals are not necessarily used but, instead, the other sensor signals such as ABS wheel speed sensor signals are used.

With the availability of the IMU sensor cluster equipped with a vehicle stability control system, several effects that are impossible to differentiate in the traditional stability control may now be included. For example, the effect of the road bank was not accurately included in the existing control algorithm. Many known systems either rely upon basic assumptions regarding conditions such as driving on a flat surface (no pitch or bank angle) or on an estimated road bank information that is usually accurate in non-event situation (e.g., steady state driving) but inaccurate in vehicle stability control events. Due to the use of INU sensor cluster in the ISS system, the road bank and grade information may now be fairly accurately determined.

Notice that the vehicle's attitudes are required to separate the gravity contamination in the acceleration sensor signals. This can be seen from the following example. The vehicle's lateral sliding used in LSC can be characterized by its lateral velocity defined along the lateral direction of the vehicle body. Such a velocity cannot be directly measured and it is determined from the lateral accelerometer measurement. The total output of the lateral accelerometer contains information that is related to the variables other than the lateral velocity. This information includes gravity contamination, the centripetal accelerations, the derivative of the lateral velocity, the sensor uncertainties. On a banked road, the gravity contributes to the lateral accelerometer measurement as significantly as the combination of the vehicle's true lateral velocity derivative and the centripetal acceleration. Due to the fact that the gravity is fixed in both its magnitude and its direction with respect to the sea level, the vehicle global attitudes can be used to find the relative position between the gravity vector and the vehicle's body directions. For this reason, the vehicle global attitudes can be used to compensate the gravity influence in the measured lateral acceleration such that the vehicle lateral velocity might be isolated and determined from the lateral acceleration sensor measurement. The same argument would hold for using pitch attitude to compensate the longitudinal accelerometer sensor signals for computing the vehicle's longitudinal velocity.

While driving on a level ground, the vehicle lateral acceleration tends to be larger than the value sustained by the limits of adhesion on a dry surface with high adhesion due to the large load transfer. On a slippery surface with low adhesion, the vehicle's lateral acceleration might be close to the limit of adhesion of the road surface. The acceleration sensors used in the existing YSC system could not provide information to differentiate this especially when a vehicle experiences large roll and pitch accelerations in excessive maneuvers. Due to the relative attitude determination from the RAD unit in the ISS system, such load transfer effects may be easily characterized. Therefore a stability control system using ISS system output would be able to achieve the same physics in the low adhesion road surface with in the high adhesion road surface, and less deviation of control performance in low adhesion road from the control performance in the high adhesion road will be experienced.

It is the objective of the current disclosure to provide an integrated sensing system and use such information for an integrated stability control system to coordinate a yaw stability control function, roll stability control function and a lateral stability control function to achieve superior control performance in comparison with the existing vehicle stability control systems.

SUMMARY

The present invention's focus is on the type of vehicle dynamic control systems that control the multiple motions of the vehicle along yaw, roll and lateral directions during the driver initiated maneuvers. The multiple motions are likely to be unstable if their controls are only conducted by an ordinary driver. For this reason, such vehicle dynamics controls are also called the vehicle stability control systems.

One theme of this invention is an integrated treatment of the potential vehicle stability controls through superior vehicle's operation state determination using an integrated sensing system. That is, an integrated stability control system (ISCS) is used to stabilize individual or multiple unstable vehicle dynamics at the same time. More specifically, integration of yaw, roll and lateral stability control systems is desirable.

The present invention provides an integrated sensing system (ISS) and uses the information from such a sensing system to conduct integrated stability control system which integrates a yaw stability control system (YSC), a roll stability control (RSC) system and a lateral stability control (LSC) system together with an anti-brakelock system (ABS), a traction control system (TCS). The integrated sensing system uses signals measured from the sensor set beyond the sensor configuration used in the traditional YSC, RSC, ABS and TCS systems. The sensors used in the ISS include a six-degree of freedom inertial measurement unit (IMU), a steering wheel angle (SWA) sensor, four ABS wheel speed sensors, a master cylinder brake pressure sensor, a brake pedal force sensor, and other sensors. The powertrain and drivetrain information and the other available information are also fed into the ISS.

Notice that the traditional ESC system integrates ABS, TCS and yaw stability management functions through individual brake pressure controls (incremental to the driver requested brake pressures). The main goal of an ESC is to augment the driver's controllability of the vehicle at the limit of the vehicle dynamics for a given road surface traction and for a given set of the driver requested driving inputs. The controllability of ESC mainly focuses on the vehicle's yaw rate response with respect to the driver's intension (determined from the inputs such as the driver's steering wheel angle input, the vehicle speed, and the like).

Although various ESC systems claim to have certain lateral stability control, its effectiveness is limited due to the fact that the sensor set used in ESC cannot accurately and robustly differentiate the effect of the road bank and the vehicle's lateral instability. If the vehicle is stable in lateral direction, the bank effect may be determined. But when the vehicle is driving close to lateral instability on a medium-sized bank road, it may not be possible to separate the vehicle's lateral siding motion from the bank influence (gravity contamination) in the lateral acceleration sensor measurement. This is especially the case when the vehicle is driven on low friction road surfaces.

The vehicle's lateral sliding usually increases the vehicle dynamically unstable tendency such that the vehicle is hard to control by ordinary drivers. Hence one of the performance requirements in vehicle stability controls is to attenuate the lateral sliding as much as possible. Notice that such a performance requirement is different from car racing, where vehicle lateral sliding is sacrificed for vehicle driving speed. One reason is that race car drivers are usually capable and experienced drivers, who can handle the vehicle well even if it is experiencing large lateral sliding. The vehicle's lateral sliding used in LSC can be characterized by its lateral velocity defined along the lateral direction of the vehicle body. Such a velocity cannot be directly measured and is determined from the lateral accelerometer measurement. The total output of the lateral accelerometer contains information that is related to the variables other than the lateral velocity. This information includes gravity contamination, the centripetal accelerations, the derivative of the lateral velocity, and the sensor uncertainties. On a banked road, the gravity contributes to the lateral accelerometer measurement as significantly as the combination of the vehicle's true lateral velocity derivative and the centripetal acceleration. Due to the fact that gravity is fixed in both magnitude and direction with respect to the sea level, the vehicle global attitudes may be used to find the relative position between the gravity vector and the vehicle's body directions. For this reason, the vehicle global attitudes can be used to compensate the gravity influence in the measured lateral acceleration such that the vehicle lateral velocity might be isolated and determined from the lateral acceleration sensor measurement. The same argument would hold for using pitch attitude to compensate the longitudinal accelerometer sensor signals for computing the vehicle's longitudinal velocity.

Sideslip angle is the inverse tangent of the ratio of the lateral velocity and the longitudinal velocity of the vehicle defined on the vehicle's body fixed coordinate system. Therefore, its accurate determination must involve the accurate determination of the longitudinal and lateral velocities of the vehicle defined on the vehicle's body-fixed coordinate system.

Of course, various other sensors such as a roll rate sensor, steering angle sensor, longitudinal acceleration sensor, and a vertical acceleration signal may be combined together or independently with the above-mentioned sensors to further refine control of the system.

In a further aspect of the invention, a method of controlling a vehicle includes determining a lateral tire force, a front lateral tire force, a rear lateral tire force, and determining a linear sideslip angle from the front lateral tire force and the rear lateral tire force. The method further includes determining a linear lateral velocity in response to the linear sideslip angle and controlling the vehicle in response to the linear sideslip angle.

In a further aspect of the invention, a method of controlling a vehicle comprises determining the outputs from an integrated sensing system, using those outputs to conduct vehicle stability controls.

In a further aspect of the invention, a method of controlling a vehicle comprises integrating YSC, RSC and LSC through the output signals calculated in the integrated sensing system.

Other advantages and features of the present invention will become apparent when viewed in light of the detailed description of the preferred embodiment when taken in conjunction with the attached drawings and appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagrammatic view of a vehicle with variable vectors and coordinate frames according to the present invention.

FIG. 2 is a block diagram of a stability system according to the present invention.

FIG. 3 is a front view of an automotive vehicle illustrating various angles according to the present invention.

FIG. 4 is a side view of an automotive vehicle illustrating various variables thereon.

FIG. 5 is a top view of a vehicle according to the present invention.

FIG. 6 is a detailed block diagrammatic view of a controller according to the present invention.

FIG. 7 is a flow chart illustrating a method of controlling the vehicle according to the present invention.

FIG. 8 is another method for determining the lateral velocity of the vehicle.

FIG. 9 is a block diagrammatic view of a longitudinal velocity estimation block

FIG. 10 is a more detailed block diagrammatic view of a method for determining longitudinal velocity.

FIG. 11 is a flow chart illustrating a method for determining the longitudinal velocity of the vehicle.

FIGS. 12A, 12B and 12C illustrate a vehicle without stability control, a vehicle with traditional stability control, and a vehicle with advanced stability control according to the present invention, respectively.

DESCRIPTION OF THE PREFERRED EMBODIMENT

In the following figures the same reference numerals will be used to identify the same components. The present invention is preferably used in conjunction with vehicle control systems including but are not limited to a yaw stability control (YSC) system, a roll stability control (RSC) system, a lateral stability control (LSC), an integrated stability control system (ISCS), or a total vehicle control system for achieving fuel economy and safety and other vehicle level performances. The system is also described with respect to an integrated sensing system (ISS) which uses a centralized motion sensor cluster such as an inertial measurement unit (IMU) and the other available but decentralized sensors. Although the centralized motion sensor such as an IMU is used, the methods described here are easily transferable to using the other discrete sensors.

Referring to FIG. 1, an automotive vehicle 10 with a vehicle dynamics control or an active safety system of the present invention is illustrated with the various forces and moments thereon during a stable vehicle dynamic condition. Vehicle 10 has front right (FR) and front left (FL) wheel/tires 12A and 12B and rear right (RR) wheel/tires 13A and rear left (RL) wheel/tires 13B, respectively. The vehicle 10 may also have a number of different types of front steering systems 14a and rear steering systems 14 b, including having each of the front and rear wheels configured with a respective controllable actuator, the front and rear wheels having a conventional type system in which both of the front wheels are controlled together and both of the rear wheels are controlled together, a system having conventional front steering and independently controllable rear steering for each of the wheels, or vice versa. Generally, the vehicle has a weight represented as Mg at the center of gravity of the vehicle, where g=9.8 m/s² and M is the total mass of the vehicle.

As mentioned above, the system may also be used with the other vehicle dynamics controls such as ride and handling control systems including active/semi-active suspension systems, anti-roll bar, or the other safety systems such as airbags or passive safety devices deployed or activated upon sensing predetermined dynamic conditions of the vehicle.

The ISS unit 16 is coupled to the ISCS unit 18. The ISS unit 16 may comprise many different sensors as will be described further below. The sensors may also be used by the ISCS unit in various determinations such as to determine a wheel lifting event such as in an imminent rollover, determine various forces including normal forces at the wheels, determine a height and position of a mass, determine the instability trend of the unstable dynamics as in unstable roll or yaw motions, determine the driver's intention, determine the feedforward control commands to drive actuators, determine feedback control commands for the desired functions, and the like. The wheel speed sensors 20 are mounted at each corner of the vehicle and generate signals corresponding to the rotational speed of each wheel which are denoted as w₁,w₂,w₃,w₄ or w_(slf),w_(srf),w_(slr),w_(srr) for left-front, right-front, left-rear and right-rear wheels, respectively. The rest of the sensors used in ISS unit 16 may include the other decentralized sensors and a centralized motion sensor such as an IMU or a RSC sensor cluster mounted directly on a rigid surface of the vehicle body such as the vehicle floor or the chassis frame. The centralized sensor cluster may not directly measure the motion variables along the vehicle body fixed longitudinal, lateral or vertical directions, it only provide motion variable information of the vehicle body at the sensor location and along the sensor's local directions.

As those skilled in the art will recognize, the frame from b₁,b₂ and b₃ is called a body frame 22, whose origin is located at the center of gravity of the car body, with the b₁ corresponding to the x axis pointing forward, b₂ corresponding to the y axis pointing off the driving side (to the left), and the b₃ corresponding to the z axis pointing upward. The angular rates of the car body are denoted about their respective axes as ω_(bx) for the roll rate, ω_(by) for the pitch rate and ω_(bz) for the yaw rate. Calculations may take place in an inertial frame 24 that may be derived from the body frame 22 as described below. The longitudinal acceleration of the vehicle body is denoted as a_(bx). The lateral acceleration of the vehicle body is denoted as a_(by). The vertical acceleration of the vehicle body is denoted as a_(bz).

As mentioned before, the centralized motion sensor might have different directions than the vehicle body-fixed directions. Such a centralized sensor cluster might be mounted on a specific location of the vehicle body. The frame system s₁s₂s₃ fixed on the sensor cluster body is called a sensor frame system. The sensor frame directions are denoted as s₁,s₂ and s₃, which are the x-y-z axes of the sensor frame. Notice that the sensor frame system is not necessary the same as the vehicle body fixed frame system b₁b₂b₃. One of the reasons is that the sensor mounting errors are possible in the production environment. For practical reasons, the centralized sensor cluster like an IMU may not be mounted on the same location which is of interest for computation purposes such as the vehicle's center of gravity or the rear axle of the vehicle body. However, it is evident that an IMU sensor measures enough information to be used to numerically translate the IMU sensor output to the motion variables at any location on the vehicle body. The six outputs of the IMU sensor along its sensor frame are denoted as ω_(sx),ω_(sy),ω_(sz),a_(sx),a_(sy) and a_(sz) for the sensor cluster's roll rate, pitch rate, yaw rate, longitudinal acceleration, lateral acceleration and vertical acceleration.

The other frame used in the following discussion includes the moving road frame, as depicted in FIG. 1. The moving road frame system r₁r₂r₃ is fixed on the driven road surface, where the r₃ axis is along the average road normal direction of a plane called a moving road plane. The moving road plane is a flat plane which serves as an average representation of the four-tire/road contact patches and it is usually a fictitious plane. If the road surface is perfectly flat, then it coincides with the moving road surface. Let the four vertical coordinates of the centers of the contact patches with respect to the inertia frame or sea level be z₀,z₁,z₂,z₃ for the front-left, front-right, rear-left and rear-right corners, then the average moving road plane experienced by the driven vehicle should have the bank angle computed as

$\begin{matrix} {{{average}\mspace{14mu} {road}\mspace{14mu} {bank}} = {\frac{1}{4}\left( {\frac{z_{0} - z_{1}}{t_{f}} + \frac{z_{2} - z_{3}}{t_{r}}} \right)}} & (1) \end{matrix}$

where t_(f) and t_(r) are the half tracks of the front and rear axles, and the inclination angle computed as

$\begin{matrix} {{{average}\mspace{14mu} {road}\mspace{14mu} {slope}} = {\frac{1}{2}\left( {\frac{z_{2} - z_{0}}{b} + \frac{z_{3} - z_{1}}{b}} \right)}} & (2) \end{matrix}$

where b is the wheel base of the vehicle.

Notice that the moving road plane is moving and yawing with the vehicle body but it is not rolling and pitching with the vehicle body. The moving road frame is the right-hand orthogonal axis system r₁r₂r₃ in which the plane containing r₁ and r₂ axes coincides with the average moving road plane. That is, r₁-axis is the projection of the longitudinal axis of the vehicle body on to the average moving road plane, r₂-axis is the projection of the lateral axis of the vehicle body on to the average moving road plane, and r₃-axis points upwards which coincides with the normal direction of the average moving road plane. Notice that the moving road frame coincides with an ISO-8855's intermediate axis system if the vehicle is driven on a level ground. On a three-dimensional road, there is no counterpart defined in IS0-8855.

Referring now to FIG. 2, the ICSC unit 18 is illustrated in further detail having an ISS unit 26 used for receiving information from a number of sensors which may include a yaw rate sensor 28, a speed sensor 20, a lateral acceleration sensor 32, a vertical acceleration sensor 33, a roll angular rate sensor 34, a hand wheel sensor 35 (steering wheel within the vehicle), a longitudinal acceleration sensor 36, a pitch rate sensor 37, steering angle (of the wheels or actuator) position sensor 38 (steered wheel angle), a suspension position (height) sensor 40. It should be noted that various combinations and sub-combinations of the sensors may be used.

The sensor cluster 16 may be disposed within a single housing 43, and includes a roll rate sensor 34 generating a roll rate signal, a pitch rate sensor 37 generating a pitch rate signal, a yaw rate sensor 38 generating a yaw rate signal, a longitudinal acceleration sensor 36 generating a longitudinal acceleration signal, a lateral acceleration sensor 32 generating a lateral acceleration signal, and a vertical acceleration sensor 33 generating a vertical acceleration.

In typical aerospace application, the centralized motion sensor like IMU is always used together with global positioning system (GPS) signals. Recently, such an approach has also made into automotive applications, see for example, “Integrating INS sensors with GPS velocity measurements for continuous estimation of vehicle sideslip and tire cornering stiffness” by D. M. Bevly, R. Sheridan and J. C. Gerdes in proceeding of the 2001 American Control Conference. In the present invention, GPS is not required.

Based upon inputs from the sensors, the ISS unit 26 may feed information to the ISCS unit 44 that further drives the actions of the available actuators. Depending on the desired sensitivity of the system and various other factors, not all the sensors may be used in a commercial embodiment. The ISCS unit 44 may control an airbag 45 or a steering actuator 46A-46D at one or more of the wheels 12A, 12B, 13A, 13B of the vehicle. Also, other vehicle subsystems such as a suspension control 48 may be used for ride, handling and stability purpose.

Roll angular rate sensor 34 and pitch rate sensor 37 may be replaced by sensors sensing the height of one or more points on the vehicle relative to the road surface to sense the roll condition or lifting of the vehicle. Sensors that may be used to achieve this include but are not limited to a radar-based proximity sensor, a laser-based proximity sensor and a sonar-based proximity sensor. The roll rate sensor 34 may also be replaced by a combination of sensors such as proximity sensors to make a roll determination.

Roll rate sensor 34 and pitch rate sensor 37 may also be replaced by sensors sensing the linear or rotational relative displacement or displacement velocity of one or more of the suspension chassis components to sense the roll condition or lifting of the vehicle. This may be in addition to or in combination with suspension distance sensor 40. The suspension distance sensor 40 may be a linear height or travel sensor and a rotary height or travel sensor.

The yaw rate sensor 28, the roll rate sensor 34, the lateral acceleration sensor 32, and the longitudinal acceleration sensor 36 may be used together to determine that single wheel or two wheels of a vehicle are lifted and the quantitative information regarding the relative roll information between the vehicle body and moving road plane such as in RSC developed at Ford. Such sensors may also be used to determine normal loading associated with wheel lift.

The roll condition such as the relative roll angle of the vehicle body with respect to the road surface or with respect to the sea level may also be established by one or more of the following translational or rotational positions, velocities or accelerations of the vehicle including the roll rate sensor 34, the yaw rate sensor 28, the lateral acceleration sensor 32, the vertical acceleration sensor 33, a vehicle longitudinal acceleration sensor 36, lateral or vertical speed sensor including a wheel-based speed sensor 20, a radar-based speed sensor, a sonar-based speed sensor, a laser-based speed sensor or an optical-based speed sensor.

The ISS unit 26 may include sensing algorithms including but not limited to reference attitude and reference directional velocity determinations, global/relative attitude determination, directional velocity determination, sensor plausibility check, sensor signal conditioning, road parameter determination, and abnormal state monitoring.

The ISS unit 26 includes various control units controlling the aforementioned sensing algorithms. More specifically, these units may include: a reference signal unit 70 (reference signal generator (RSG)), which includes an attitude reference computation and a velocity reference computation, a road profile unit 72 (road profile determination unit (RPD)), an attitude unit or relative attitude determination unit 74 (RAD), a global attitude unit 76 (global attitude determination unit (GAD) and a directional velocity unit 78 (directional velocity determination unit (DVD)), a sensor plausibility unit 80 (sensor plausibility check unit (SPC)), an abnormal state unit 82 (abnormal state monitoring unit (ASM)), a sensor signal compensating unit 84 (SSC), an estimation unit 86 (force and torque estimation unit (FATE)), a car body to fixed reference frame unit 88 (body to reference unit (B2R)), a normal loading unit 90 (normal loading determination unit (NLD)), a vehicle parameter unit 92 (vehicle parameter determination unit (VPD)), a four wheel driver reference model 94 and a sideslip angle computation 96. Signals generated from any one of the aforementioned units are referred to prediction of vehicle operation states signals.

The ISCS unit 44 may control the position of the front right wheel actuator 46A, the front left wheel actuator 46B, the rear left wheel actuator 46C, and the right rear wheel actuator 46D. Although as described above, two or more of the actuators may be simultaneously controlled. For example, in a rack-and-pinion system, the two wheels coupled thereto are simultaneously controlled. Based on the inputs from sensors 28 through 43B, controller 26 determines a roll condition and/or wheel lift and controls the steering position and/or braking of the wheels.

The ISCS unit 44 may be coupled to a brake controller 60. Brake controller 60 controls the amount of brake torque at a front right brake 62 a, front left brake 62 b, rear left brake 62 c and a rear right brake 62 d. The functions performed through ISCS might include an RSC function 110, a YSC function 66, and an LSC function 69. The other functional units such as an anti-lock-braking system (ABS) unit 64 and a traction control system (TCS) unit 65 may be provided by the brake suppliers and are usually controlled through signals calculated from the supplier's sensing and control module (which might uses the same or part of the sensors used in this invention) Those functions might be improved through utilizing the signals calculated in ISS unit.

Speed sensor 20 may be one of a variety of speed sensors known to those skilled in the art. For example, a suitable speed sensor may include a sensor at every wheel that is averaged by the ISS unit 26. The algorithms used in ISS may translate the wheel speeds into the travel speed of the vehicle, Yaw rate, steering angle, wheel speed, and possibly a slip angle estimate at each wheel may be translated back to the speed of the vehicle at the center of gravity. Various other algorithms are known to those skilled in the art. Speed may also be obtained from a transmission sensor. For example, if speed is determined while speeding up or braking around a corner, the lowest or highest wheel speed may not be used because of its error. Also, a transmission sensor may be used to determine vehicle speed instead of using wheel speed sensors.

Although the above discussions are valid for general stability controls, some specific considerations of using them in RSC will be discussed. The roll condition of a vehicle during an imminent rollover may be characterized by the relative roll angle between the vehicle body and the wheel axle and the wheel departure angle (between the wheel axle and the average road surface). Both the relative roll angle and the wheel departure angle may be calculated in relative roll angle estimation module (RAD 74) by using the roll rate, lateral acceleration sensor signals and the other available sensor signals used in ISS unit. If both the relative roll angle and the wheel departure angles are large enough, the vehicle may be in either single wheel lifting or double wheel lifting. On the other hand, if the magnitude of both angles is small enough, the wheels are likely all grounded, therefore the vehicle is not rolling over. In case that both of them are not small and the double wheel lifting condition is detected or determined, the sum of those two angles will be used to compute the feedback commands for the desired actuators so as to achieve roll(over prevention. The variables used for this purpose might be included in the ISS unit.

The roll information of a vehicle during an imminent rollover may be characterized by rolling radius-based wheel departure roll angle, which captures the angle between the wheel axle and the average road surface through the dynamic rolling radii of the left and right wheels when both of the wheels are grounded. Since the computation of the rolling radius is related to the wheel speed and the linear velocity of the wheel, such rolling-radius based wheel departure angle will assume abnormal values when there are large wheel slips. This happens when a wheel is lifted and there is torque applied to the wheel. Therefore, if this rolling radius-based wheel departure angle is increasing rapidly, the vehicle might have lifted wheels. Small magnitude of this angle indicates the wheels are all grounded. The variables used for this purpose might be included in the ISS unit.

The roll condition of the vehicle during an imminent rollover may be seen indirectly from the wheel longitudinal slip. If during a normal braking or driving torque the wheels at one side of the vehicle experience increased magnitude of slip, then the wheels of that side are losing longitudinal road torque. This implies that the wheels are either driven on a low mu surface or lifted up. The low mu surface condition and wheel-lifted-up condition may be further differentiated based on the chassis roll angle computation, i.e., in low mu surface, the chassis roll angle is usually very small. The variables used for this purpose might be included in the ISS unit.

The roll condition of the vehicle during an imminent rollover may be characterized by the normal loading sustained at each wheel. Theoretically, when a normal loading at a wheel decreases to zero, the wheel is no longer contacting the road surface. In this case a potential rollover is underway. Large magnitude of this loading indicates that the wheel is grounded. Normal loading is a function of the calculated chassis roll and pitch angles. The variables used for this purpose might be included in the ISS unit.

The roll condition of a vehicle during imminent rollover may be identified by checking the actual road torques applied to the wheels and the road torques, which are needed to sustain the wheels when they are grounded. The actual road torques may be obtained through torque balancing for each wheel using wheel acceleration, driving torque and braking torque. If the wheel is contacting the road surface, the calculated actual road torques must match or be larger than the torques determined from the nonlinear torques calculated from the normal loading and the longitudinal slip at each wheel. The variables used for this purpose might be included in the ISS unit.

The roll condition of a vehicle during an imminent rollover may be characterized by the chassis roll angle itself, i.e., the relative roll angle between the vehicle body and the wheel axle. If this chassis roll angle is increasing rapidly, the vehicle might be on the edge of wheel lifting or rollover. Small magnitude of this angle indicates the wheels are not lifted or are all grounded. Therefore, an accurate determination of the chassis roll angle is beneficial for determining if the vehicle is in non-rollover events and such computation is conducted in RAD unit in ISS.

The roll condition of a vehicle during imminent rollover may also be characterized by the roll angle between the wheel axle and the average road surface, which is called a wheel departure angle (WDA). If the roll angle is increasing rapidly, the vehicle has lifted wheel or wheels and aggressive control action needs to be taken in order to prevent the vehicle from rolling over. Small magnitude of this angle indicates the wheels are not lifted. The variables used for this purpose might be included in the ISS unit.

The ISCS unit 18 may include control function logic and control function priority logic. As illustrated, the logic resides within ISCS unit 44, but may be part of the ISS unit 26 and/or brake controller 60.

Referring now to FIG. 3, the relationship of the various angles of the vehicle 10 (end view of which is shown) relative to the road surface 11 is illustrated. In the following, a reference road bank angle θ_(bank) is shown relative to the vehicle 10 on a road surface. The vehicle has a vehicle body 10 a and wheel axle 10 b. The wheel departure angle θ_(wda) is the angle between the wheel axle and the road. The relative roll angle θ_(xr) is the angle between the wheel axle 10 b and the body 10 a. The global roll angle θ_(x) is the angle between the horizontal plane (e.g., at sea level) and the vehicle body 10 a.

Another angle of importance is the linear bank angle. The linear bank angle is a bank angle that is calculated more frequently (perhaps in every loop) by subtracting the relative roll angle generated from a linear roll dynamics of a vehicle (see U.S. Pat. No. 6,556,908 which is incorporated by reference herein), from the calculated global roll angle (as one in U.S. Pat. No. 6,631,317 which is incorporated by reference herein). If all things were slowly changing without drifts, errors or the like, the linear bank angle and reference road bank angle terms would be equivalent.

During an event causing the vehicle to roll, the vehicle body is subject to a roll moment due to the coupling of the lateral tire force and the lateral acceleration applied to the center of gravity of vehicle body. This roll moment causes suspension height variation, which in turn results in a vehicle relative roll angle (also called chassis roll angle or suspension roll angle). The relative roll angle is used as an input to the activation criteria and to construct the feedback brake pressure command for RSC function, since it captures the relative roll between the vehicle body and the axle. The sum of such a chassis roll angle and the roll angle between wheel axle and the road surface (called wheel departure angle) provides the roll angle between the vehicle body and the average road surface, which is feeding back to the RSC module.

Referring now to FIGS. 4 and 5, an automotive vehicle 10 is illustrated with various parameters illustrated thereon. The side view of automotive vehicle 10 is illustrated. A front suspension 82 and a rear suspension 82 r are illustrated. The suspensions are coupled to the body at a respective suspension point 84 f, 84 r. The distance from the suspension point 84 f to the center of the wheel is labeled z_(sh). The distance from the center of gravity CO to the front suspension is labeled as b_(f). The distance from the CG to the rear suspension point 84 r is labeled as b_(r). The vertical distance between the center of gravity and the suspension point are labeled respectively as h_(f) and h_(r). A portion of the body axis b₃ and the road axis r₃ are illustrated. The angle therebetween is the relative pitch angle θ_(yr), The rolling radius of the tire is labeled as z_(w).

Referring specifically now to FIG. 5, a top view of vehicle 10. Lateral and longitudinal velocities of the center of gravity are denoted as v_(x) and v_(y), a yaw angular rate is denoted as ω_(z), a front wheel steering angle is denoted as δ_(s), lateral acceleration is represented by a_(y), and longitudinal acceleration is represented by a_(x). Also illustrated is the front track width t_(f) and rear track width t_(r). That is, the track widths are half the track widths from the center line or center of gravity of the vehicle to the center of the tire. Rear axle 14b extends between the rear wheels 13 a, 13 b. Various calculations may be performed relative to the rear axles of 14 b.

Referring now to FIG. 6, the interrelationship among the various units in a controlled vehicle dynamics is shown. The vehicle dynamics is represented by the module 102. The vehicle dynamics includes the angular and translation movements of the vehicle. The motion and actuation sensors 28-40 generate signals corresponding to the vehicle dynamics and the actions of the various actuators. The sensors are fed into the ISS unit 104. The sensor signals and the calculated signals from a system other than ISS (for example, the brake supplier's own brake control computations) 106 may be fed into to sensing module 108. The ISS unit is coupled to the ISCS unit 44 and may specifically be coupled to but not limited to the RSC function module 110 and a YSC function module 66. A function priority system command 80 may also be included. The output of the TCS module 68 and the ABS module 64 may be coupled to priority logic system command 114, which in turn is coupled to the actuation distribution and commanding block 116. The actuation distribution commanding block 116 is coupled to the powertrain control ECU 118 and to the brakes 119. The powertrain control module 118 may be coupled to the engine 120 and transmission 122. The actuation of the engine, the brakes and the transmission may affect the vehicle dynamics 102, which in turn is sensed by the various sensors. Thus, as can be seen, a continuous loop of sensing and controlling with respect to the vehicle dynamics is illustrated.

Referring now to FIG. 7, further detail of the operation ISS unit 26 is illustrated. The various computations in the ISS unit begin in start step 300. In step 302, signals from the various sensors are received. The sensor electronics usually provide an indication of the quality of the sensor signal sent out from a specific sensor and the environmental temperature around the sensors may also be provided. Built-in self-test of the sensors may also be conducted through the sensor electronics. In step 302, the signals from various sensors and sensor electronics are received. In this step, some extra actions other than “receiving” may be conducted. For example, the sensor signals are monitored for electric failures that may include communication failures, wiring failures and power supply failures. The communication failures may be due to a defective CAN hardware device, high frequency interference, CAN bus off events, missing CAN message from the sensor electronics, too high CAN load, wrong CAN controller configuration, and the like. The wiring failures might be due to interruptions, short circuit, electronics interference and leakage. The power supply failures might be due to short circuit to ground or to the wrong voltage, low voltage of the power supply. The sensor signals are also monitored for sudden changes and for out-of range monitoring. A sudden signal change with very high change rate might be a good indication of electronic failure of the sensor.

In step 304, the plausibility of the sensors is determined in the sensor plausibility check module (SPC). The main focus of SPC is on sensor defects, mechanical failures, or an implausible operation content. For example, a loose sensor might have certain implausible sensor (vibration) content which does not match the vehicle's motion; a sensor might be mounted incorrectly on the vehicle; the sensor electronics is gradually degraded (e.g., electronics component aging, electronic wiring aging) that cannot be detected through electronic monitoring). In all these cases, the overall vehicle dynamics behavior is checked to determine if certain sensor signals are plausible. For example, a roll rate plausibility check is conducted by comparing the roll rate sensor value with an roll rate estimation based on sensor signals other than roll rate such as the lateral acceleration, the yaw rate, the wheel speed and the steering wheel signals; a yaw rate plausibility check is conducted by comparing the yaw rate sensor value with an estimated yaw rate value based on the steering wheel angle, the lateral acceleration and wheel speed signals.

In step 30C, if at least one sensor is not plausible, step 308 is executed in which sensor fault processing logic is performed. The plausibility flags are set and saved in an accessible memory. A reduced model computation will be conducted and a partial function control will be determined. The info from the sensor fault processing logic will be fed to the ISCS unit 44 through the signal flow 309.

In step 310, sensor signals are compensated in the sensor signal compensation unit SSC. SSC will conduct sensor zero-point offset compensation, sensor signal drift (for example, due to ambient temperature change) compensation), the sensor noise compensation, etc. The sensor mounting error such as sensor misaligiments (see U.S. Pat. No. 6,782,315) might also be computed at this module.

In step 312, the abnormal operation states of the vehicle subsystems are monitored in the abnormal state monitoring module ASM. The abnormal states of the vehicle might be a flat tire or a tire with significant inflation pressure drop, a significant loading bias (such as a full trunk loading in a truck or a large roof loading in a SUV), off-shelf wheels, extreme wear of the brake pads, the imbalanced wheel (see U.S. Pat. No. 6,834,222). In step 313, if there are abnormal states presented, the abnormal state processing unit is conducted in step 314. The abnormal state processing logic involves assessing the severity of the abnormal states in terms of influence of the ISCS functions. If the abnormal operation states are not critical, then certain compensations of the involved signals will be compensated and the control function in ISCS will be continued. For example, if ASM detects one of the tires has a large tire pressure drop (but still above the threshold for a flat tire), the wheel's rolling radius will be learned and the wheel speed of that wheel will be compensated for by the learned rolling radius. The same is true for large size tires, spare tire, roof loading and trunk loading. If the abnormal states are critical such as in a flat tire case, the normal vehicle function like ISCS will be partially or completely shut off.

If there are no abnormal states detected in step 313, the system proceeds to step 316 to compute the reference signals which will be used in reference signal generator RSG (see U.S. Pat. No. 7,010,409). Notice that the variables calculated in RSG may be used by various functional units such as SPC, GAD, DVD, etc.

In step 318, the vehicle body's roll and pitch attitudes with respect to the sea level are calculated in GAD module. Although the engineering definition of the vehicle body's roll and pitch attitudes are pretty obvious (which are simultaneously measured just like in the previous definitions of the road bank and slope), there are no mathematical characterizations existed. Mathematically, Euler used sequential rotation to define the angular distance between different frame systems. If the longitudinal axis is used for the roll axis, and lateral axis for the pitch axis, then the Euler angles along those two axes are called the roll and pitch Euler angles respectively. Experimentally, the engineering definition (simultaneous measurement) of the vehicle body's global roll and pitch angles are very close to the Euler attitudes (sequential rotational measurement) and the differences are negligible. For this reason, the roll and pitch angles defined here are interchangeable with the Euler angles.

The Euler angles of the sensor cluster with respect to the sea level can be related to the roll, pitch and yaw angular rate sensor signals through the following coupled kinematics

{dot over (θ)}_(sx)=ω_(sx)+[ω_(sy) sin(θ_(sx))+ω_(sz) cos(θ_(sx))]tan(θ_(sy))

{dot over (θ)}_(sy)=ω_(sy) cos(θ_(sx))−ω_(sz) sin(θ_(sx))   (3)

The relationship depicted above reveals complicated nonlinear functions and it indicates that a simple integration of the roll rate could provide accurate information about the roll attitude only if both the pitch and yaw rate are negligible, which means the vehicle is dominated by roll motion; the roll attitude angle and yaw rate are negligible, the pitch attitude is limited; and, the pitch attitude angle is negligible with non-extreme pitch and yaw rates.

Similarly, the simple integration of the pitch rate could lead to accurate prediction of the pitch attitude angle if the roll attitude angle is negligible and the yaw rate is not extremely large. If the acceleration projected from the sensor frame to the vehicle body frame is defined as a_(bx), a_(by), the vehicle directional velocities measured on the vehicle body frame as v_(bx), v_(by), and the attitudes of the vehicle body with respect to the sea level as θ_(bx), θ_(by), then the following equations 4 hold.

{dot over (v)} _(bx) =a _(bx)+ω_(bz) v _(by) +g sin θ_(by)

{dot over (v)} _(by) =a _(by)−ω_(bz) v _(bx) −g sin θ_(bx) cos θ_(by)   (4)

In step 320, the vehicle body's relative attitudes with respect to its axles or the moving road frame are calculated. Those attitudes might be computed as in U.S. Pat. Nos. 6,556,908, 6,684,140, and 6,631,317.

In step 322, the signals measured along the sensor frames are transformed to signals along the moving road frames in the B2R module. If an IMU sensor cluster is used, the sensor signals might need to be compensated for frame different. For example, the sensor signals are meant to be along the sensor frame directions are not necessary along the vehicle body frame directions or the directions of control interests such as the moving road frames. The angular differences between different frames can be mathematically characterized using Euler transformations and using Euler angles. The Euler angles are defined as the sequential rotation along an axis to move one frame system to another frame system. When the Euler angle between the sensor frame s₁s₂s₃ and the body frame b₁b₂b₃ is Δθ_(x),Δθ_(y),Δθ_(z), then the angular rates and accelerations of the vehicle body defined along the body-fixed frame b₁b₂b₃ but at the sensor location may be calculated using the following matrix computation

$\begin{matrix} {\begin{bmatrix} \omega_{bx} \\ \omega_{by} \\ \omega_{bz} \end{bmatrix} = {{{{E\left( {{\Delta\theta}_{z},{\Delta\theta}_{y},{\Delta\theta}_{x}} \right)}\begin{bmatrix} \omega_{sx} \\ \omega_{sy} \\ \omega_{sz} \end{bmatrix}}\begin{bmatrix} a_{bx} \\ a_{by} \\ a_{bz} \end{bmatrix}} = {{E\left( {{\Delta\theta}_{z},{\Delta\theta}_{y},{\Delta\theta}_{x}} \right)}\begin{bmatrix} a_{sx} \\ a_{sy} \\ a_{sz} \end{bmatrix}}}} & (5) \end{matrix}$

where the matrix E(,,) is defined as the following for three Euler angles θ₃,θ₂,θ₁

${E\left( {\theta_{3},\theta_{2},\theta_{1}} \right)} = {\quad\begin{bmatrix} {\cos \; \theta_{3}\cos \; \theta_{2}} & {{\cos \; \theta_{3}\sin \; \theta_{2}\sin \; \theta_{1}} - {\sin \; \theta_{3}\cos \; \theta_{1}}} & {{\cos \; \theta_{3}\sin \; \theta_{2}\cos \; \theta_{1}} + {\sin \; \theta_{3}\sin \; \theta_{1}}} \\ {\sin \; \theta_{3}\cos \; \theta_{2}} & {{\cos \; \theta_{3}\cos \; \theta_{1}} + {\sin \; \theta_{3}\sin \; \theta_{2}\sin \; \theta_{1}}} & {{\sin \; \theta_{3}\sin \; \theta_{2}\cos \; \theta_{1}} - {\cos \; \theta_{3}\sin \; \theta_{1}}} \\ {{- \sin}\; \theta_{2}} & {\cos \; \theta_{2}\sin \; \theta_{1}} & {\cos \; \theta_{2}\cos \; \theta_{1}} \end{bmatrix}}$

Notice that the difference between the sensor frame and the body frame is mainly due to the sensor mounting errors. The sensor cluster mounting direction errors may be obtained by converting the simultaneously measured distances at the fixed points among the two frame systems to angular distance as in the definition of the moving road surface before. Mathematically, such a simultaneous measurement is different from the sequential definition of the Euler angles. However the differences between those two measurements are not significant. For this reason, the Δθ_(x),Δθ_(y),Δθ_(z) used in this invention is interchangeable with the engineering terminology such as the roll misalignment, pitch misalignment and yaw misalignment. U.S. Pat. No. 6,782,315 provides a method for determining Δθ_(x),Δθ_(y),Δθ_(z). If the vehicle body motion signals projected along the moving road frame is of interest, a further Euler transformation can be conducted. When the engineering roll angle between the vehicle body and the axle is θ_(xy) (which can be defined by using the suspension height differences between the left and right sides of the wheels and which is the same as the aforementioned relative or chassis roll angle, see U.S. Pat. Nos. 6,556,908 and 6,684,140), the engineering pitch angle between the vehicle body and the axle is θ_(yr) (which can be defined by using the suspension height differences between the front and rear axles and which is the same as the aforementioned relative or chassis pitch angle, see U.S. Pat. Nos. 6,556,908 and 6,684,140), the engineering roll angle between the axle and the moving road frame is θ_(wda) (which is non-zero only if the vehicle is experiencing an imminent rollover, see US patent publications 2004/0064236 and 2004/0162654), then the vehicle body's motion variables ω_(mrfx),ω_(mrfy),ω_(mrfz),a_(mrfx),a_(mrfy),a_(mrfz) along the moving road frames may be computed as in the following

$\begin{matrix} {\begin{bmatrix} \omega_{mrfx} \\ \omega_{mrfy} \\ \omega_{mrfz} \end{bmatrix} = {{{{E\left( {0,\theta_{yr},{\theta_{xr} + \theta_{wda}}} \right)}\begin{bmatrix} \omega_{bx} \\ \omega_{by} \\ \omega_{bz} \end{bmatrix}}\begin{bmatrix} a_{mrfx} \\ a_{mrfy} \\ a_{mrfz} \end{bmatrix}} = {{E\left( {0,\theta_{yr},{\theta_{xr} + \theta_{wda}}} \right)}\begin{bmatrix} a_{bx} \\ a_{by} \\ a_{bz} \end{bmatrix}}}} & (6) \end{matrix}$

By using the features of the Euler transformation, ω_(mrfx),ω_(mrfy),ω_(mrfz),a_(mrfx),a_(mrfy),a_(mrfz) can be directly related to the sensor signals as in the following

$\begin{matrix} {\begin{bmatrix} \omega_{mrfx} \\ \omega_{mrfy} \\ \omega_{mrfz} \end{bmatrix} = {{{{E\left( {{\Delta\theta}_{z},{{\Delta\theta}_{y} + \theta_{yr}},{{\Delta\theta}_{x} + \theta_{xr} + \theta_{wda}}} \right)}\begin{bmatrix} \omega_{bx} \\ \omega_{by} \\ \omega_{bz} \end{bmatrix}}\begin{bmatrix} a_{mrfx} \\ a_{mrfy} \\ a_{mrfz} \end{bmatrix}} = {{E\left( {{\Delta\theta}_{z},{{\Delta\theta}_{y} + \theta_{yr}},{{\Delta\theta}_{x} + \theta_{xr} + \theta_{wda}}} \right)}\begin{bmatrix} a_{bx} \\ a_{by} \\ a_{bz} \end{bmatrix}}}} & (7) \end{matrix}$

Notice that in many cases, the motion variables such as accelerations at locations other than the sensor location are of interests. For example, in this invention disclosure, the vehicle body accelerations a_(mrfrax),a_(mrfray),a_(mrfraz) at the vehicle rear axle location are of interests for lateral stability control and they can be further transformed according to the following computations:

a _(mrfrax) =a _(mrfx) +p _(xsl2ra)(ω_(mrfy) ²+ωmrfz²)

a _(mrfray) =a _(mrfx) −p _(xsl2ra)({dot over (ω)}_(mrfz)+ω_(mrfx)ω_(mrfy))

a _(mrfraz) =a _(mrfz) +p _(xsl2ra)({dot over (ω)}_(mrfy) −ω _(mrfx)ω_(mrfz))   (8)

where p_(xsl2ra) is the distance between the IMU sensor location and the rear axle in the longitudinal direction. The sensor is assumed to be located on the center line of the vehicle body. If such an assumption is not true, extra terms will need to be added to the above computation.

In 324, the vehicle's directional velocities such as the vehicle longitudinal velocity, lateral velocity, are calculated. The detailed description about such computations will be provided below.

In step 326, the vehicle's sideslip angle is computed. The detailed description about such a computation will be provided below.

In 328, the lateral and longitudinal tire forces, braking torques and driving torques applied to the wheels are computed in the FATE unit by using the available motion sensor information such as body motion variables from IMU sensors, wheel speed sensors and actuator specific information such as the caliper brake pressure and the engine axle torque.

In 330, the normal loading from the road applied to the wheels at the wheel contact patches are computed in the NLD unit.

In 332, the vehicle's operational parameters such as the vehicle roll gradient, the vehicle loading and mass, the vehicle's c.g. height, etc., are determined in the VPD unit.

In 334, the driving road parameters such as the road bank, slope, friction level, etc., are determined in the RPD unit. Some of those computations might be found in U.S. Pat. No. 6,718,248 and US patent publication 2004/0030475.

In 336, the driver's driving references are determined based on the driver reference model DRM. One of the DRMs might be a bicycle model used in finding the vehicle yaw reference, which is used to determine the driver's intension. Notice that the tire forces reside on the moving road plane, hence a driving reference model using tire force balancing must work everything in the moving road frame. Let the desired motion variables of the vehicle body at the vehicle body center of gravity location but projected to the moving road frame a_(mrfcgxd),a_(mrfcgyd),v_(mrfcgyd),ω_(mrfcgzd) for longitudinal acceleration, lateral acceleration, lateral velocity and yaw rate, satisfy the following relationship based on a four-wheel vehicle model (here only the front wheel steering case is presented, for rear wheel steering and four wheel steering case, similar equations may be obtained)

I _(z){dot over (ω)}_(mrfcgzd) =F _(x1) l cos(γ+δ_(s))−F _(x2) l cos(γ−δ_(s))+(F _(x3) −F _(x4))t ₄ +F _(y1) l sin(γ+δ_(s))+F _(y2) l sin(γ−δ_(s))−F _(yr) l _(r)

M ₁ a _(mrfcgxd)=(F _(x1) +F _(x2))cos(δ_(s))+F _(x3) +F _(x4)−(F _(y1) +F _(y2))sin(δ_(s))

M ₁ a _(mrfcgyd)=(F _(x1) +F _(x2))sin(δ_(s))+(F _(y1) +F _(y2))cos(δ_(s))+F _(yr)   (9)

where l=√{square root over (l_(f) ²+l_(f) ²)}, γ=a tan(t_(f) ll_(f)). M₁ is the total vehicle mass; a_(mrfcgxd) and a_(mrfcgyd) are the desired vehicle c.g. longitudinal and lateral acceleration projected on the moving road frame; δ_(s) is the wheel steering angle, which may be calculated through the driver's steering wheel angle and the known steering gear ratio; F_(y1), F_(y2) and F_(yr) are the lateral forces applied to the front-left, front-right, rear axle (sum of the lateral forces of the rear-left and the rear right wheels) from the road. F_(x1),F_(x2),F_(x3),F_(x4) are the four longitudinal tire forces at the front-left, the front-right, the rear-left and the rear-right wheels. If the desired front wheels' sideslip angle is α_(mrffd) and the rear wheels' sideslip angle is α_(mrfrd) (both of them are defined on the moving road plane) then the desired vehicle sideslip angle, the steering inputs and the desired yaw rate are related with each other as in the following:

$\begin{matrix} {{\alpha_{mrffd} = {\frac{v_{mrfcgyd}}{v_{mrfcgx}} + \frac{l_{f}\omega_{mrfcgzd}}{v_{mrfcgx}} - \delta_{S}}}{\alpha_{mrfrd} = {\frac{v_{mrfcgyd}}{v_{mrfcgx}} - \frac{l_{r}\omega_{mrfcgzd}}{v_{mrfcgx}}}}} & (10) \end{matrix}$

Using those wheel sideslip angles, the lateral tire forces may be expressed as

F_(y1)=c₁α_(mrffd)

F_(y2)=c₂α_(mrffd)

F_(yr)=c_(r)α_(mrfrd)   (11)

where c₁ is the cornering stiffness of the vehicle and by using those forces, then

I _(z){dot over (ω)}_(mrfcgzd) =a(δ_(s) , v _(mrfcgx))ω_(mrfcgzd) +b(δ_(s) , v _(mrfcgx))v _(mrfcgyd) +c(δ_(s))F _(x)

M ₁ {dot over (v)} _(mrfcgyd) =e(δ_(s) , v _(mrfcgx))ω_(mrfcgzd) +f(δ_(s) , v _(mrfcgx))v _(mrfcgyd) +g(δ_(s))F _(x) +h(θ_(mrfx))   (12)

The above defined desired yaw rate and the desired lateral velocity needs to satisfy the following equation for longitudinal direction force balancing

y(F _(x),δ_(y) ,a _(mrfcgx))=c _(f) x(v _(mrfcgyd),ω_(mrfcgzd) v _(mrfcgx),δ_(s))   (13)

where c_(f)=c₁+c₂, and functions y() and x() may be expressed as the following

$\begin{matrix} {{x = {\left( {\frac{v_{mrfcgyd}}{v_{mrfcgx}} + \frac{l_{f}\omega_{mrfcgzd}}{v_{mrfcgx}} - \delta_{S}} \right){\sin \left( \delta_{S} \right)}}}{y = {{\left( {F_{x\; 1} + F_{x\; 2}} \right){\cos \left( \delta_{S} \right)}} + F_{x\; 3} + F_{x\; 4} - {M_{f}a_{mrfcgx}}}}} & (14) \end{matrix}$

and this equation may be used to adaptively estimate the average cornering stiffness experienced by the vehicle under the desired motion as in the following iterative scheme

c _(f)(k+1)=c _(f)(k)−ρw(k+1)[y(k)−c _(f)(k)x(k)]ΔT   (15)

where ρ is the adaptive gain and ΔT is the sampling time of the computation. Notice that in the above DRM computation, the inputs are the vehicle's longitudinal velocity projected on the moving road frame but computed at the c.g. of the vehicle, the measured longitudinal acceleration of the vehicle body projected on the moving road frame at c.g. of the vehicle body, the driver's steering wheel inputs. Those variables are directly related to the sensor measurements. The outputs are the desired yaw rate ω_(mrfcgzd), the desired lateral velocity v_(mrfcgyd), the updated tire cornering stiffness c_(f).

All the above calculated signals are fed into the ICSC module 44 and the arbitrated or prioritized or overridden control commands are then sent to the ECU or hardware electronics of the available actuators through the control signal represented by 342. The specific applications using those signals will be discussed in the following. The ISCS according to a first embodiment of the present invention intends to use differential braking to keep the vehicle moving in stable conditions regardless the driver's inputs or the road surface conditions. The instability of the vehicle during a driver-initiated maneuver could be in roll direction (such as a RSC event), yaw direction (such as a YSC event), lateral direction (such as a LSC event), or any combination of RSC, YSC and LSC events.

For the LSC events case, the lateral instability of the vehicle happens when its lateral excursion relative to the driver's desired path is divergent. For turning purpose, the vehicle has to have certain (absolute) lateral excursion (especially in the front axle) which is defined through the driver's steering input. Even if the driver is requesting a large steering input (which means a large absolute lateral excursion), the vehicle could still be laterally stable (since the vehicle's relative lateral excursion could be small). During lateral instability, the velocity of this relative lateral excursion is non-zero and without crossing zero for a significant amount of time (for example, for more than 1 or 2 seconds). Mathematically, when a vehicle is negotiating a path, the total vehicle's velocity coincides with the tangent direction of the path and the vehicle is laterally stable. On the other hand, if the vehicle's total velocity largely deviates from the tangent direction of the path defined by the driver's steering input, the vehicle is likely in the laterally instable mode. Considering the rear wheel does not have any steering angle (for a front wheel steering vehicle), when the center of the rear axle negotiates the path, its total velocity must be along the longitudinal direction of the moving road frame. For this reason, the angle between the total velocity of the center of the rear axle and the vehicle's longitudinal direction are used as the measurements for the vehicle's lateral instability. The rear location is used in this invention as the location to define the vehicle's lateral instability. Since the vehicle's path always resides on the road surface or the moving road frame, the lateral instability also needs to use the angle between the total velocity of the rear axle projected on the moving road frame and the longitudinal direction of the moving road frame. Such an angle is a variation of the normally used sideslip angle at the rear and is set forth as the variable β_(LSC). β_(LSC) as defined as in the following:

$\begin{matrix} {\beta_{LSC} = {\beta_{mrfra} = {\tan^{- 1}\left( \frac{v_{mrfray}}{v_{mrfrax}} \right)}}} & (16) \end{matrix}$

Ideally, the desired value for β_(LSC) is zero. For practical implementation purpose, a non-zero desired value β_(LSCd) is used. The error signal

β_(LSCe)=β_(LSC)−β_(LSCd)

is used as a feedback signal to compute the control command for conducting LSC control. Let β_(db1) and β_(db2) be the two dead-bands for compute the torques necessary for counteracting the vehicle's lateral instability. The proportional portion of the feedback torque M_(LSCP) can be computed as in the following

if during oversteer condition  if β_(LSCe) ≧ β_(db1) and β_(LSCe) ≦ β_(db2)   M_(LSCP) = 0  else if β_(LSCe) < β_(db1) (18)   M_(LSCP) = G_(SP)(β_(db1) − β_(LSCe))  else   M_(LSCP) = G_(SP)(β_(db2) − β_(LSCe)) end where G_(sp) is a control gain which could be an adaptive gain based on the vehicle's dynamics. The over-steer condition of the vehicle can be determined based on the steering wheel information, the lateral acceleration information, the yaw rate information, the yaw error information, the sideslip angle information, etc. One example of such an over-steer condition can be determined from the YSC yaw rate error. The variable dβ_(LSCe) is the derivative of β_(LSCe), which can be directly computed from β_(LSCe) or from the other signals. The variables dβ_(db1) and β_(db2) are the two dead-bands. The derivative portion of the feedback torque M_(LSCD) can be computed as in the following:

if during oversteer condition  if d β_(LSCe) ≧ d β_(db1) and d β_(LSCe) ≦ d β_(db2)   M_(LSCD) = 0  else if d β_(LSCe) < d β_(db1) (19)   M_(LSCD) = G_(SD)(d β_(db1) − d β_(LSCe))  else   M_(LSCD) = G_(SD)(d β_(LSCe) − d β_(LSCe)) end where G_(SD) is a control gain which could be an adaptive gain based on the vehicle's dynamics. The under-steer condition of the vehicle can be determined based on the steering wheel information, the lateral acceleration information, the yaw rate information, the yaw error information, the sideslip angle information and the like. One example of an under-steer condition may be determined from the YSC yaw rate error. The final feedback torque M_(LSC) for LSC control can be obtained through the following blending

M _(LSC)=λ₁ M _(LSCP)+λ₂ M _(LSCD)   (20)

where the weights λ₁,λ₂ are two positive numbers less than 1 which weight the relative importance of the proportional term and the derivative term in the final feedback torque. In certain driving conditions, a constraint may be set forth as λ₁=λ₂=1.

The YSC yaw moment demand M_(YSC) may be generated through yaw rate error feedback control. When ω_(mrfcgzd) is the desired or target yaw rate which is determined in DRM module, ω_(mrfz) is the vehicle body's yaw rate but projected on the moving road frame. Then the vehicle yaw rate error can be calculated as

ω_(mrfze)=ω_(mrfz)−ω_(mrfcgzd)   (21)

If ω_(mrfze)>ω_(zosthrehold)>0 together with the sign of the vehicle's sideslip angle matching the turning direction of the vehicle, where ω_(zosthrehold) is a variable that reflects the threshold of the vehicle's over-steer characteristics, the vehicle is likely to be in an over-steer condition, i.e., the vehicle's yaw response exceeds the driver's intention determined from the driver's steering wheel angle. Based on this, the following counteracting torque applied to the yaw direction of the vehicle may be computed so as to reduce the discrepancy between the actual yaw rate and the desired yaw rate

if (ω_(mrfe) − ω_(zosthreshold) ≧ ω_(os) & &a_(mrfray) ≧ a_(yos) & & β_(mrfra) ≦ −β_(os)  || ω_(zosthreshold) − ω_(mrfe) ≧ ω_(os) & &a_(mrfray) ≦ −a_(yos) 0 & & β_(mrfra) ≧ β_(os))   M_(YSC) = G_(OSYP)(ω_(zosthreshold) − ω_(mrfe)) + (22)   G_(OSYD)(dω_(zosthreshold) − ω_(mrfe)); else   M_(YSC) = 0 where ω_(os),a_(yos),β_(os) are the dead-bands for the yaw rate error beyond threshold ω_(zosthreshold), the dead-band for the lateral acceleration, the dead-band for the sideslip angle at the rear axle; dω_(zosthreshold),{dot over (ω)}_(mrfe) are the derivatives of the threshold ω_(zosthreshold) and the yaw rate error; G_(YPOS),G_(YDOS) are feedback control gains which are functions of the driving conditions. M_(OS) is the control torque along the yaw direction that tries to regulate the yaw error during over-steer condition.

If ω_(mrfze)<ω_(zosthrehold)<0 with ω_(zosthrehold) together with the sign of the vehicle's sideslip angle matching the turning direction of the vehicle, where ω_(zosthrehold) is a variable that reflects the threshold of the vehicle's under-steer characteristics, the vehicle is likely to be in a under-steer characteristics, i.e., the vehicle's yaw response is below the driver's intension determined from the driver's steering wheel angle

f (ω_(mrfe) − ω_(zosthreshold) ≦ −ω_(us) & &a_(mrfray) ≧ a_(yus)  || ω_(zosthreshold) − ω_(mrfe) ≦ ω_(us) & &a_(mrfray) ≦ −a_(yus))   M _(US) = G_(YPUS)(ω_(zosthreshold) − ω_(mrfe)) + (23)   G_(YDUS)(dω_(zosthreshold) − ω_(mrfe)); else   M _(US) = 0 where ω_(us),a_(yus),β_(us) are the dead-bands for the yaw rate error beyond threshold ω_(zosthreshold), the dead-band for the lateral acceleration, the dead-band for the sideslip angle at the rear axle; dω_(zosthreshold) is the derivative of the threshold ω_(zosthreshold); G_(YPUS),G_(YDUS) are feedback control gains which are functions of the driving conditions. M_(US) is the control torque along the yaw direction that tries to regulate the yaw error during under-steer condition.

The RSC control is an over-steer type of control. The RSC control command can be computed from the yaw torque based on the feedback signal of vehicle body-to-road roll angle. This torque may be set forth as M_(RSC).

Since for the different functions, there are different computations for commanding the control torques applied to the vehicle's yaw direction through differential braking. Those control torques regulate different error signals such as the sideslip angle at the rear axle of the vehicle but projected on the moving road plane (for LSC), the yaw error of the vehicle along the normal direction of the moving road plane (for YSC) and the relative roll angle between the vehicle body and the moving plane (for RSC) There is a need to arbitrate or prioritize or integrate those different computations to achieve the optimized integrated stability control performance. For under-steer control, since there is only one control function (in YSC under-steer control), such an integration is not necessary. However, for over-steer control, such an integration is required. If M_(OS) and M_(US) are defined as the finally arbitrated over-steer control torque and the under-steer control torque, then one of such integration scheme can be summarized as in the following, where the weights γ₁,γ₂ and γ₃ are three positive numbers less than 1 that weight the relative importance among the three feedback torques:

if (M_(RSC) > 0 & &M_(LSC) > 0 & & M_(YSC) > 0 || M_(RSC) < 0 & &M_(LSC) < 0 & & M_(YSC) < 0) {  M_(OS) = γ₁M_(RSC) + γ₂M_(LSC) + γ₃M_(YSC);M_(US) = 0; } elseif (| M_(RSC) |> p_(rscth) & & | M_(LSC) |< p_(lscth) & & | M_(YSC) |< p_(yscth)) {  M_(OS) = M_(RSC);M_(US) = 0; } elseif (| M_(RSC) |< p_(rscth) & & | M_(LSC) |> p_(lscth) & & | M_(YSC) |< p_(yscth)) {  M_(OS) = M_(LSC);M_(US) = 0; } elseif (| M_(RSC) |< p_(rscth) & & | M_(LSC) |< p_(lscth) & & | M_(YSC) |> p_(yscth)) {  M_(OS) = M_(YSC);M_(US) = 0; } elseif (| M_(RSC) |> p_(rscth) & & | M_(LSC) |> p_(lscth) & & | M_(YSC) |< p_(yscth)) {  [M_(OS) ,i] = max(| M_(YSC) |,| M_(LSC) |);M_(OS) = M_(OS)sign(M_(i));M_(US) = 0; } elseif (| M_(RSC) |> p_(rscth) & & | M_(LSC) |< p_(lscth) & & | M_(YSC) |> p_(yscth)) {  [M_(OS) ,i] = max(| M_(RSC) |,| M_(YSC) |);M_(OS) = M_(OS)sign(M_(i));M_(US) = 0; } elseif (| M_(RSC) |< p_(rscth) & & | M_(LSC) |> p_(lscth) & & | M_(YSC) |> p_(yscth)) {  [M_(OS) ,i] = max(| M_(LSC) |,| M_(YSC) |);M_(OS) = M_(OS)sign(M_(i));M_(US) = 0; } else { (24)  [M_(OS) ,i] = max(| M_(RSC) |,| M_(LSC) |,|M_(YSC) |);M_(OS) = M_(OS)sign(M_(i));M_(US) = 0; }

After the afore-mentioned torque M_(OS) or M_(US) has been computed, it may be realized by building or releasing braking pressure on selective wheels or wheel. Generally speaking, in over-steer case, the controlled wheel is the front outside wheel or both the front outside wheel and the rear outside wheel, i.e., the brake pressure modification is requested for the front outside wheel or for the outside wheels; in under-steer control, the controlled wheel is the rear inside wheel or both the inside wheels, i.e., brake pressure modification is requested for the rear inside wheel or for the inside wheels. Based on the aforementioned rule of thumb, various types of controls can be performed based on the various driving conditions. In the following discussion, a rear wheel driving vehicle is assumed and the similar consideration holds for front wheel driving or four wheel driving vehicles.

In the case of under-steering with throttle open, the engine torque is first reduced so as to regulate the under-steer yaw error. If the engine torque reduction is not enough to bring the vehicle's yaw error below certain threshold, the rear inside wheel is braked and the brake pressure amount is computed based on the yaw error feedback torque M_(US). For example, the brake pressure for the inside rear wheel needs to generate the following longitudinal slip ratio

$\begin{matrix} {\rho_{ir} = \frac{M_{US}}{t_{r}\kappa_{ir}N_{ir}}} & (25) \end{matrix}$

where ρ_(ir) is the longitudinal slip-ratio of the inside rear wheel, t_(r) is the half track of the rear axle, κ_(ir),N_(ir) are the gain and the normal loading for the wheel at the inside rear location. If the requested insider rear wheel slip ratio |ρ_(ir)|≧ ρ with ρ being the allowed minimum slip ratio at the rear wheels to prevent the rear end to loose. In some cases, when the vehicle's sideslip angle β_(LSC) is kept below certain threshold, both front inside and rear inside wheels might be applied brake pressure based on the following requested slip ratios

$\begin{matrix} {{{\rho_{ir} = \overset{\_}{\rho}};}{\rho_{if} = \frac{M_{US} - {\overset{\_}{\rho}t_{r}\kappa_{if}N_{if}}}{t_{r}\kappa_{if}N_{if}}}} & (26) \end{matrix}$

If after the above brake based under-steer control, the vehicle's under-steer yaw error is still big and the driver's steering input is still increasing, the vehicle might have significant tire lateral force saturation in the front wheels. In this case, the vehicle may be slowed down through additional braking. Such additional braking could be applied to all the four wheels. Since the vehicle's turning radius is proportional to the vehicle's speed, hence slowing down the vehicle is indirectly helping reduce the under-steer of the vehicle. Notice that whenever the vehicle sideslip angle is beyond certain threshold, the afore-mentioned under-steer control request will be zeroed out. If the turning-off is not quick enough to change the actual caliper pressure, both under-steer control and over-steer control might happen briefly together. If the brake pressure apply to two wheels at the same axle, M_(US) is might be realized through longitudinal force differential as in the following, the yaw moment demand is developed through the side-to-side deviation of the longitudinal slip ratio between the two wheels

$\begin{matrix} {{\Delta\rho} = \frac{M_{US} - {t_{f}\left( {{\kappa_{1}N_{1}\lambda_{fl}} - {\kappa_{2}N_{2}\lambda_{fl}}} \right)}}{t_{f}\kappa_{2}N_{2}}} & (27) \end{matrix}$

In the case of under-steering with vehicle coasting, then the brake pressure will be first applied to the rear inside wheel. If the under-steer yaw rate error cannot be significantly reduced, then a further increase of the brake pressure in the front inside wheel (as in 26)) will be conducted. If after the above brake based under-steer control, the vehicle's under-steer yaw error is still big enough, further braking will be applied to the outside wheels so as to further slow down the vehicle. Notice that whenever the vehicle sideslip angle is beyond certain threshold, the afore-mentioned under-steer control request will be zeroed out.

In the case of under-steering with driver braking, the braking pressure is sent to the front inside wheel and the rear inside wheel until ABS is activated. If the rear inside wheel enters ABS and the vehicle is still demanding large yaw error correction command, then part of the pressure for the rear inside wheel will be sent to the front inside wheel. If the front inside wheel enters ABS and the rear inside wheel is not in ABS and the yaw command is still high, part of the front inside wheel pressure will be redirected to the rear inside wheel. Notice that whenever the vehicle sideslip angle is beyond certain threshold, the afore-mentioned under-steer control request will be zeroed out. If such a turning-off is not quick enough to change the actual caliper pressure, both under-steer control and over-steer control might happen briefly together.

In the case of over-steering with throttle open, the engine torque reduction is also conducted and at the same time a brake pressure based on M_(OS) is sent to the front outside wheel.

In the case of over-steering with vehicle coasting, a brake pressure based on M_(OS) is sent to the front outside wheel.

In the case of over-steering with driver braking, the braking pressure is applied to the front outside wheel and the front outside wheel is allowed to lock so as to kill significant yaw error and sideslip angle.

In the case of over-steering with driver braking, applying pressure to the front outside wheel. If the front outside wheel is locked but the yaw command is still large, the rear inside wheel slip control is activated (reducing the rear inside wheel slip ratio so as to increase the lateral force at rear axle).

More detailed control scheme to realized the calculated feedback torque M_(OS) and M_(US) can be similarly developed based on the drivetrain configuration, etc.

The further details of the computation used in ISS will be provided here.

The kinematics used in the ISS unit can be described by the following equations for the body frame

dθ _(bx)=ω_(bx)+[ω_(by) sin θ_(bx)+ω_(bz) cos θ_(bx)]tan θ_(by)

dθ _(by)=ω_(by) cos θ_(bx)−ω_(bz) sin θ_(by)

dv _(bx) =a _(bx)+ω_(bz) v _(by) +g sin θ_(by)

dv _(bray) =a _(bray)−ω_(bz) v _(bx) −g sin θ_(bx) cos θ_(by)   (28)

where v_(bray),a_(bray) are the lateral velocity and acceleration at the rear axle location but projected on vehicle body-fixed frames. dθ_(bx),dθ_(by),dv_(bx),dv_(bray) are the variables whose theoretic values are the time derivatives of θ_(bx),θ_(by),v_(bx),v_(bray), but in reality due to sensor error or sensor uncertainties, they are not necessarily equal to {dot over (θ)}_(bx),{dot over (θ)}_(by),{dot over (v)}_(bx),{dot over (v)}_(bray). In the sensor frame, the above equations can be rewritten as

dθ _(sx)=ω_(sx)+[ω_(sy) sin θ_(sx)+ω_(sz) cos θ_(sx)]tan θ_(sy)

dθ _(sy)=ω_(sy) cos θ_(sx)−ω_(sz) sin θ_(sx)

dv _(sx) =a _(sx)+ω_(sz) v _(sy) +g sin θ_(sy)

dv _(sray) =a _(sray)−ω_(sz) v _(sx) −g sin θ_(sx) cos θ_(sy)   (29)

In the moving road frame, the following equations are true

dθ _(mrfx)=ω_(mrfx)+[ω_(mrfy) sin θ_(mrfx)+ω_(mrfz) cos θ_(mrfx)]tan θ_(mry)

dθ _(mrfy)=ω_(mrfy) cos θ_(mrfx)−ω_(mrfz) sin θ_(mrfx)

dv _(mrfx) −a _(mrfx)+ω_(mrfz) v _(mrfy) +g sin θ_(mrfy)

dv _(mrfray) =a _(mrfray)−ω_(mrfz) v _(mrfx) −g sin θ_(mrfx) cos θ_(mrfy)   (30)

In the following discussion, the details of directional velocity determination in ISS will be provided. First, start from the computation of the linear lateral velocity based on the bicycle model is performed.

Linear Lateral Velocity

Referring now to FIG. 8, the traditional bicycle model may be used to determine, the yaw and lateral dynamics of a vehicle. This may be performed using the following differential equations

I _(z){dot over (ω)}_(mrfz) =b _(f) F _(mrff) cos δ−b _(r) F _(mrfr)−(I _(y) −I _(x))ω_(mrfy)ω_(mrfx) +Y _(z)

M ₁ a _(mrfcgy) =F _(mrff) cos δ+F _(mrfr)   (31)

where F_(mrff) and F_(mrfr) are the lateral forces applied to the front axle and rear axle through tires respectively; Y_(z) is the yawing moment due to yaw stability control, which may be estimated based on the desired yaw stability command and the estimated road surface μ; I_(x),I_(y),I_(z) are the roll, pitch and yaw moments of inertia of the vehicle along the axes of the moving road frame which may be approximated by those in the vehicle body frame; M, is the vehicle total mass; b_(f) is the distance from the vehicle center of gravity to the front axle; b_(r) is the distance from the vehicle center of gravity to the rear axle; δ is the steering angle at the front wheels.

Let β_(mrff) and β_(mrfr) be the front and rear wheel sideslip angles respectively, which may be expressed as

$\begin{matrix} {{\beta_{mrflinf} = {\beta_{mrflincg} + \frac{b_{f}\omega_{mrfz}}{v_{mrfx}} - \delta}}{\beta_{mrflinr} = {\beta_{mrflincg} - \frac{b_{r}\omega_{mrfz}}{v_{mrfx}}}}} & (32) \end{matrix}$

For small wheel sideslip angles, i.e.

|β_(mrff)|≦ β _(f)

|β_(mrfr)|≦ β _(r)   (33)

the lateral forces applied to the front axle and the rear axle may be approximated as the following linear relationship

F_(mrff)≈−c_(f)β_(mrff)

F_(mrfr)≈−c_(r)β_(mrfr)   (34)

where β _(f) and β _(r) are front and rear wheel slip angle upper bounds to guarantee the linear relationships in (34), and c_(f) and c_(r) are the cornering stiffness of the front and rear tires respectively. c_(f) and c_(r) usually vary with the road condition and the vehicle operation condition. Due to the steering angle and large loading towards the front axle, the front wheel sideslip angle β_(mrff) usually exceeds the linear range of the tire force much earlier than the rear axle tires. Notice as soon as (33) is satisfied, (31)-(34) are true. If the vehicle is operated in its nonlinear dynamics range, using (34) with fixed values for c_(j) and c_(r) is no longer close to the true wheel sideslip angles. However the computation may still be used for other computation purpose, for example, used in the switch logic for other computations. For this reason, the following linear sideslip angles may be defined regardless of whether the vehicle is operated in its linear dynamics range or its nonlinear dynamics range

$\begin{matrix} {{\beta_{mrflinf} = \frac{F_{mrff}}{c_{f}}}{\beta_{mrflinr} = \frac{F_{mrfr}}{c_{r}}}} & (35) \end{matrix}$

where the lateral tire forces may be calculated as in the following, which are true for all the driving conditions (independent of road surface mu level, the aggressiveness of the driver's driving input, and if the tires entering their nonlinear ranges, etc.)

$\begin{matrix} {{F_{mrff} = \frac{{I_{z}{\overset{.}{\omega}}_{mrfz}} + {\left( {I_{y} - I_{x}} \right)\omega_{mrfy}\omega_{mrfx}} - Y_{z} + {b_{r}M_{t}a_{mrfy}}}{\left( {b_{f} + b_{r}} \right)\cos \; \delta}}{F_{mrfr} = \frac{{{- I_{z}}{\overset{.}{\omega}}_{mrfz}} - {\left( {I_{y} - I_{x}} \right)\omega_{mrfy}\omega_{mrfx}} + Y_{z} + {b_{f}M_{t}a_{mrfy}}}{\left( {b_{f} + b_{r}} \right)}}} & (36) \end{matrix}$

In step 700 the sensor information may be transformed into a body frame of reference. The raw rear wheel sideslip angle β_(mrfinrraw) is determined in step 702 and is used here as an estimation and control reference. It may be further expressed as the following:

β_(mrflinrraw) =p _(yawgrad) a _(mrfy) −pyawaccoef{dot over (ω)}_(mrfz) −p _(coupling)ω_(mrfy)ω_(mrfx) +p _(yawtorq) Y _(z)   (37)

where the four composite parameters, called yaw gradient, yaw acceleration coefficient, coupling coefficient and yaw torque gain respectively, may be related to the vehicle parameters as in the following

$\begin{matrix} {{P_{yawfrad} = \frac{b_{f}M_{t}}{\left( {b_{f} + b_{r}} \right)c_{r}}}{p_{yawacccoef} = \frac{I_{z}}{\left( {b_{f} + b_{r}} \right)c_{r}}}{p_{coupling} = \frac{\left( {I_{y} - I_{x}} \right)}{\left( {{b\;}_{f} + b_{r}} \right)c_{r}}}{p_{yawtorq} = \frac{1}{\left( {b_{f} + b_{r}} \right)c_{r}}}} & (38) \end{matrix}$

The four composite parameters may also be estimated through the fitting of the empirical vehicle data to (37) during linear vehicle dynamics. This can also be done in real-time adaptively.

Those parameters may also be calculated based on vehicle testing data. For example, using a Least Square Parameter Identification for the maneuvers within the linear tire force range, the calculated β_(mrflinrraw) may be used to fit the measured sideslip angle.

In order to take account of the effect lateral load transfer, a load transfer correction factor is delivered in step 704.

if (|θ_(xr) |≦0.75 p_(rollgrad)) { ${\beta_{mrflinrraw} = \frac{\beta_{mrflinrraw}}{1 - {p_{csxmax}\; \left( {{\theta_{xr}}/p_{rollgrad}} \right)^{2}}}};$ } (39) else { ${\beta_{mrflinrraw} = \frac{\beta_{mrflinrraw}}{1 - {0.5625\mspace{11mu} p_{csxmax}}}};$ }

The final linear sideslip angle may be calculated in step 706 in the following:

β_(mrflinr) =p _(yawmdl)β_(mrflmr)+(1−p _(yawmdl))β_(mrflinrraw)   (40)

where θ_(xr) is the relative roll angle calculated in RSC, which reflects the lateral load transfer, p_(rollgrad) is the roll gradient and p_(csxmax) is the maximum percentage change of the cornering stiffness reduction due to lateral load transfer. The linear lateral velocity based on the linear sideslip angle in step 708 is defined as

v _(mrflinray) =v _(mrfx) tan(β_(mrflinr))   (41)

If an equivalent cornering stiffness is defined at the rear wheel c_(eqr) in the following

$\begin{matrix} {c_{eqr} = \frac{F_{mrfr}}{\beta_{mrfr}}} & (42) \end{matrix}$

the equivalent cornering stiffness ratio in step 710 may be computed as

$\begin{matrix} {S = \frac{c_{r}}{c_{eqr}}} & (43) \end{matrix}$

Since c_(r) is proportional to the surface friction level μ (where μ=1) and c_(eqr) is proportional to the actual road surface μ level, the equivalent cornering stiffness ratio indirectly relates to the road surface μ, although the actual function would be rather complicated.

If the vehicle's true (nonlinear) sideslip angle β_(mrfr) is known, then the aforementioned equivalent cornering stiffness ratio may be calculated as in the following:

$\begin{matrix} {S = \frac{\beta_{mrfr}}{\beta_{mrflinr}}} & (44) \end{matrix}$

However, β_(mrfr) is not known and needs to be calculated first, but {dot over (β)}_(mrfr) is known.

A good computation strategy would be working on {dot over (β)}_(mrfr) and using washout filter to obtain some rough computation of S, and using such a calculated S in some other conditions related to the road surface μ. The raw value S_(raw) of the equivalent cornering stiffness ratio S is computed through a washout filter and the computation may be summarized as in the following

if (|β_(mrflinr)|>p_(minlatvel)) v_(temp)= β_(mrflinr)v_(mrfx); else if (β_(mrflinr)>=0 && β_(mrflinr) ≦ p_(minlatvel)) v_(temp)= p_(minlatvel)v_(mrfx); (45) else v_(temp)=−p_(minlatvel)v_(mrfx); ${S_{raw} = {{\left\lbrack {\gamma_{sintcoef} - \frac{\left( {{{\overset{.}{v}}_{mrfx}\beta_{mrflinr}} + {v_{mrfx}{\overset{.}{\beta}}_{mrflinr}}} \right){\Delta T}}{v_{temp}}} \right\rbrack S_{raw}} + \frac{{\overset{.}{v}}_{mrfy}{\Delta T}}{v_{temp}}}};$ where γ_(sintcoef) is the washout filter coefficient, p_(minlatvel) is the lower bound of the linear sideslip angle β_(mrflinr), ΔT is the sampling time, v_(mrfx), is the longitudinal velocity of the vehicle.

Considering at low vehicle speeds, the above computation may be contaminated by the signal noises and the other errors, hence at low vehicle speed, the equivalent cornering stiffness ratio is set to 1. Since S also reflects the road surface μ, thus during known road surface cases, S needs to be limited by the known value. Therefore,

if (v_(mrfx) < p_(smspd) ) S_(raw) = 1; else { if ( μ_(max) > p_(smmu)) (46) ${S_{raw} = {\min \; \left( {\frac{1}{p_{smmu}},{\max \; \left( {S_{raw},1} \right)}} \right)}};$ else ${S_{raw} = {\min \; \left( {\frac{1}{p_{smmu}},{\max \; \left( {S_{raw},1} \right)}} \right)}};$ }

The final equivalent cornering stiffness ratio is obtained by feeding the computed raw value through a low-pass filter in the following:

S=d _(lpf1) S+n _(lpf1) S _(raw)   (47)

where n_(lpf1),d_(lpf1) are coefficients for the low pass filter.

Lateral Velocity and Its Reference

DVD schemes are set forth below. The first, second, or combination of both, may be used. The first scheme is illustrated in step 711 a and uses a longitudinal kinematics constraint. If the longitudinal velocity is known, then the third equation in (29) or (30) becomes redundant. That is, the third equation in (29) or (30) may be eliminated. However, a total elimination of such an equation would not be desirable since potential useful information may be lost. A case where the vehicle is driven steady state on a constant radius turn on a level ground is illustrated. In such a driving case, the acceleration balancing through the first velocity equation in (29) or (30) may actually be used to constrain the lateral velocity (the so-called centrifugal acceleration)

$\begin{matrix} {{a_{mrfx} = {{- v_{mrfy}}\omega_{mrfz}}}{or}{v_{mrfy} = {- \frac{a_{mrfx}}{\omega_{mrfz}}}}} & (48) \end{matrix}$

Therefore, the first velocity equation may be used as a constraint for the estimated lateral velocity. Such a constraint is called the longitudinal kinematics constraint, which is expressed as

k _(x)(v _(mrfray) ,v _(mrfraz))=ā _(mrfx)ω_(mrfz) v _(mrfray)−ω_(mrfy) v _(mrfraz) {dot over (v)} _(mrfx)   (49)

Notice that the true lateral velocity and longitudinal velocity should satisfy

k _(x)(v _(mrfray) ,v _(mrfraz))=0   (50)

and good estimations {circumflex over (v)}_(mrfx) and {circumflex over (v)}_(mrfray) of v_(mrfx) and v_(mrfray) should make |k_(x)(v_(mrfray),v_(mrfraz))| as small as possible The estimation error is:

e _(mrfy) =v _(mrfray) −{circumflex over (v)} _(mrfray)

e _(mrfz) =v _(mrfraz) −{circumflex over (v)} _(mrfraz)   (51)

By using the kinematics constraint shown in (51), the following estimation scheme may be formulated as:

{circumflex over ({dot over (v)} _(mrfray) −ā _(mrfray)−ω_(mrfz) v _(mrfx)+ω_(mrfx) {circumflex over (v)} _(mrfx) +G ₁ k _(x)({circumflex over (v)} _(mrfray) ,{circumflex over (v)} _(mrfraz))

{dot over ({circumflex over (v)} _(mrfraz) =ā _(mrfraz)ω_(mrfy) v _(mrfx)−ω_(mrfx) {circumflex over (v)} _(mrfray) +G ₂ {circumflex over (k)} _(x)({circumflex over (v)} _(mrfray) ,{circumflex over (v)} _(mrfraz))   (52)

where G₁ and G₂ are gains needed to be properly chosen. The error dynamics for estimation scheme (52) satisfies

ė _(mrfy)=ω_(mrfx) e _(mrfz) −G ₁(ω_(mrfy) e _(mrfz)−ω_(mrfz) e _(mrfy))

ė _(mrfz)=−ω_(mrfx) e _(mrfy) −G ₂(ω_(mrfy) e _(mrfz)−ω_(mrfz) e _(mrfy))   (53)

The gains G₁ and G₂ may then be selected such that the errors in (59) approach zero. If the following function is chosen

$\begin{matrix} {V = {\frac{1}{2}\left( {e_{mrfy}^{2} + e_{mrfz}^{2}} \right)}} & (54) \end{matrix}$

then its derivative along the error dynamics (53) may be calculated as

{dot over (V)}=G ₁ω_(mrfz) e _(mrfy) ² −G ₂ω_(mrfz) −G ₁ω_(mrfy) e _(mrfy) +G ₂ω_(mrfz) e _(mrfy) e _(mrfz)   (55)

If a positive number ρ and the gains G₁ and G₂ are chosen in the following

G ₁=−ρω_(mrfz)

G₂=ρω_(mrfy)   (56)

then (55) can be made negative

{dot over (V)}=−ρ(ω _(mrfz) e _(mrfy)−ω_(mrfy) e _(mrfz))²≦0   (57)

which implies the choice of gains in (56) will keep the rate change of the estimation error decreasing. Hence, (52) provides a useful construction of the lateral velocity, which utilizes feedback information.

If the lateral angle Θ_(λσ) is defined as:

$\begin{matrix} {\Theta_{x\; \sigma} = {\int_{t_{1}}^{\sigma}{{\omega_{mrfx}(\tau)}\ {\tau}}}} & (58) \end{matrix}$

and considering v_(mrfraz) is close to zero and is negligible, (52) may be further simplified

$\begin{matrix} {{v_{mrfray}\left( t_{2} \right)} = {\sec \; \Theta_{{xt}_{2}}\left\{ {{v_{mrfray}\left( t_{1} \right)} + {\int_{t_{1}}^{t_{2}}{\left\lbrack {{\left( {{\overset{\_}{a}}_{mrfray} - {\omega_{mrfz}v_{mrfx}} - {\rho \; \omega_{sz}{k_{x}\left( {{\hat{v}}_{mrfray},{\hat{v}}_{mrfraz}} \right)}}} \right)\cos \; \Theta_{x\; \sigma}} - {\left( {{\overset{\_}{a}}_{mrfraz} + {\omega_{mrfy}v_{mrfx}} + {\rho \; \omega_{sy}{k_{x}\left( {{\hat{v}}_{mrfray},{\hat{v}}_{mrfraz}} \right)}}} \right)\sin \mspace{7mu} \Theta_{x\; \sigma}}} \right\rbrack {\sigma}}}} \right\}}} & (59) \end{matrix}$

In the second scheme set forth, a road constraint alignment (RCA) in step 711 b is determined. An RCA presumes the vehicle cannot lift off from the ground for a long period of time. Hence, the vertical velocity must be very small and its derivative must not have significant low frequency content. If the vertical velocity is neglected, then the vertical velocity equation becomes redundant. However, a total elimination of such an equation would be bad since potential useful information is lost. Considering a case where the vehicle is in a dynamic turn on a smooth road surface and the vehicle has no-zero roll motion, the vehicle lateral velocity may be constrained in the following equation (the so-called centrifugal acceleration):

$\begin{matrix} {{{\omega_{mrfx}v_{mrfray}} = {a_{mrfraz} + {\omega_{mrfy}v_{mrfx}}}}{or}{v_{mrfy} = \frac{a_{mrfraz} + {\omega_{mrfy}v_{mrfx}}}{\omega_{mrfx}}}} & (60) \end{matrix}$

Generally speaking, RCA can be described as:

k _(z)(v _(mrfx) ,v _(mrfray))=LPF[ā _(mrfraz)ω_(mrfy) v _(mrfx)−ω_(mrfx) v _(mrfray)]=0   (61)

where LPF stands for the low-pass filtering. This may be referred to as a vertical kinematics constraint. For estimated velocities, (61) may be thought of as an enforcement mechanism for compensating potential error sources. Denote the estimated longitudinal and vertical velocities as {circumflex over (v)}_(mrfx) and {circumflex over (v)}_(mrfray), the estimation errors as

e _(mrfx) =v _(mrfx) −{circumflex over (v)} _(mrfx)

e _(mrfy) =v _(mrfray) −{circumflex over (v)} _(mrfray)   (62)

One of the estimation schemes may be described as in the following:

{dot over ({circumflex over (v)}_(mrfx) =ā _(mrfx)+ω_(mrfz) {circumflex over (v)} _(mrfray) +H ₁({circumflex over (v)} _(mrfx) −v _(mrfx))+H ₂ k _(z)({circumflex over (v)} _(mrfx) ,{circumflex over (v)} _(mrfray))

{dot over ({circumflex over (v)} _(mrfray) =ā _(mrfray)−ω_(mrfz) {circumflex over (v)} _(mrfx) +H ₃ k _(z)({circumflex over (v)} _(mrfx) ,{circumflex over (v)} _(mrfray))   (63)

where H₁, H₂ and H₃ are gains needed to be properly chosen. Notice that v_(sx) is calculated from wheel speed sensor signals together with certain calculated signals, which means it is independent of computations in (65). The error dynamics for the estimation scheme (63) satisfies

ė _(mrfx)=ω_(mrfz) e _(mrfx) −H ₁ e _(mrfx) −H ₂ LPF[ω _(mrfx) e _(mrfy)−ω_(mrfy) e _(mrfx)]

ė _(mrfy)=−ω_(mrfz) e _(mrfx) −H ₂ LPF[ω _(mrfx) e _(mrfy)−ω_(mrfy) e _(mrfx)]  (64)

For analysis simplicity, the LPF may be removed.

Let

$\begin{matrix} {V = {\frac{1}{2}\left( {e_{mrfx}^{2} + e_{mrfy}^{2}} \right)}} & (65) \end{matrix}$

then the derivative of V along the error dynamics (65) may be calculated as

{dot over (V)}−H ₁ e _(mrfx) ² −H ₂ω_(mrfx) e _(mrfy) e _(mrfx) +H ₂ω_(mrfy) e _(mrfx) ² −H ₃ω_(mrfx) e _(mrfy) ² +H ₃ω_(mrfy) e _(mrfx)   (66)

If

H₁>0

H ₂=−ρω_(mrfy)

H₃=ρω_(mrfx)   (67)

then

{dot over (V)}=H ₁ e _(mrfx) ²−ρ(ω_(mrfz) e _(mrfy)−ω_(mrfy) e _(mrfz))²≦0   (68)

which implies that the choice of gains in (67) will keep rate change of the estimation error decrease By replacing the estimated longitudinal velocity with the external longitudinal velocity (for example, from wheel speed sensors), the following provides a construction of the intended estimation of the lateral velocity

$\begin{matrix} {{{v_{mrfray}\left( t_{2} \right)} = {{{- {v_{mrfx}\left( t_{2} \right)}}\; \tan \; \left( \Theta_{{zt}_{2}} \right)} + {{\sec\left( \; \Theta_{{zt}_{2}} \right)}\left\{ {{v_{mrfray}\left( t_{1} \right)} + {\int_{t_{1}}^{t_{2}}{\left\lbrack {{\left( {{\overset{\_}{a}}_{mrfx} - {\rho \; \omega_{mrfy}{k_{z}\left( {{\hat{v}}_{mrfx},{\hat{v}}_{mrfray}} \right)}}} \right){\sin\left( \; \Theta_{z\; \sigma} \right)}} + {\left( {{\overset{\_}{a}}_{mrfray} + {\rho \; \omega_{mrfz}{k_{z}\left( {{\hat{v}}_{mrfx},{\hat{v}}_{mrfray}} \right)}}} \right)\cos \; \left( \ \Theta_{z\; \sigma} \right)}} \right\rbrack {\sigma}}}} \right\}}}}{where}} & (69) \\ {\Theta_{z\; t} = {\int_{t_{1}}^{t}{{\omega_{mrfz}(\tau)}\ {\tau}}}} & (70) \end{matrix}$

Notice that the above two schemes can be used only when the dynamics measured by the IMU sensors experiences large angular rates. If the angular rates are small, the afore-mentioned constraints is no longer useful. For this reason, the computation described before cannot be directly used to compute the lateral velocity of the vehicle body at the rear axle but projected on the moving road frame. However, they could serve as a reference lateral velocity.

In step 712, the reference lateral velocity is computed. Considering vertical velocity is small, the vertical velocity equation in the previous RCA can be dropped, Therefore,

dv _(mrfray) =ā _(mrfray)−ω_(mrfz) v _(mrfx)−ρω_(mrfz)(ā _(mrfx) +ω _(mrfz) {circumflex over (v)} _(mrfray) −{dot over (v)} _(mrfx))   (71)

Equation (71) is accurate enough when the vehicle does not have significant heave and roll motion. If the vehicle is driven on a bumpy road with significant roll motion, the cross product term ω_(mrfx)v_(mrfraz) needs to be estimated and add to (71) so as to remove the error.

Equation (71) may be integrated to obtain a signal which is a reference signal for the lateral velocity at the rear axle and projected on the moving road frame. Due to the computational errors in the attitudes and the sensor errors, both the integration rate and the initial condition will be changed based on certain reliable computations.

One of such reliable signals is the linear sideslip angle of the vehicle at its rear axle, which is already calculated above. The magnitude of such a variable may be calculated a in the following

β_(mrflinrmag) =p _(dipflbetam)β_(mrflinrmag)+(1−p _(dipflbetam)|β_(mrflinr)|)   (72)

The roll attitude of the moving road frame θ_(mrfx) may have computational errors, the integration coefficient used will be adapted based on the magnitude of θ_(mrfx) as in the following

if (| θ_(mrfx) |≦ θ₁)  γ_(intcoef) = p_(intcoef1); else if (θ₁ ≦| θ_(mrfx) |≦ θ₂)  γ_(intcoef) = p_(intcoef2); else if (θ₂ ≦| θ_(mrfx) |≦ θ₃) (73)  γ_(intcoef) = p_(intcoef3); else if (θ₂ ≦| θ_(mrfx) |≦ θ₃)  γ_(intcoef) = p_(intcoef4); else  γ_(intcoef) = p_(intcoef5);

The further adaptation of the integration rate is conducted based on the indication of the road surface μ reflected by the maximum acceleration of the vehicle and the magnitude of the equivalent cornering stiffness ratio computed above

$\mu_{\max} = {{n_{{lpf}\; 2}\mu_{\max}} + {\frac{d_{{lpf}\; 2}}{g}\sqrt{a_{mrfrax}^{2} + a_{mrfray}^{2} + a_{mrfraz}^{2}}}}$ if (  (μ_(max) > p_(muminhighmu) or S > p_(smaxhighmu) ) and  |a_(mrfray)|≦ p_(minlata) ) (74) { γ_(intcoef) = min(γ_(intcoef), γ_(intcoeffast)); }

The conditional integration of (71) may be conducted in the following

if ( ( ( β_(mrflinrmag) ≧ p_(smlinb) && S ≦ p_(minslowmu) )  or ω_(mag) ≧ p_(smrat) || | d_(swa) |> p_(smswa)) & &v_(mrfx) ≧ p_(smspd) ) { k_(cgain) = k_(gain)(v_(mrfx) ≧ p_(gainswitchspd) );  k_(cx) = a_(mrfrax) + w_(mrfz)v_(mrfrefy) − v_(mrfx);  if(k_(cx) ≧ k_(cxmax) ) (75)   k_(cx) = k_(cxmax);  if(k_(cx) ≦ −k_(cxmax) )   k_(cx) = −k_(cxmax);  v_(mrfrefy) = v_(mrfrefy)γ_(intcoef) + (a_(mrfray) − w_(mrfz)v_(mrfx) −  k_(cgain)w_(mrfz)k_(cx)); }

Let the final computation of the lateral velocity at the rear axle be v_(mrfray). The above calculated lateral velocity reference may be bounded by the final computation of v_(mrfray). The following conditioning will be conducted if the final computation of v_(mrfray) exceeds the magnitude of the final computation v_(mrfrefy), the reference value will be gradually adjusted to approach the final computation v_(mrfray)

if ( ( v_(mrfrefy) ≧ v_(mrfray) and v_(mrfray) ≧ 0 )  or ( v_(mrfrefy) ≧ v_(mrfray) and v_(mrfray) ≦ 0 ) ) (76)  v_(mrfrefy) = v_(mrfray) + (v_(mrfrefy) − v_(mrfray))/ p_(rho);

A further conditioning will be conducted at linear sideslip range, which assigns the computation to the linear velocity calculated from the linear sideslip angle

if ( ( v_(mrfrefy) v_(mrflinray) < 0 or (v_(mrflinr) ≧ 0 and v_(mrfrefy) ≦ v_(mrflinray))  or (v_(mrflinr) ≦ 0 and v_(mrfrefy) ≧ v_(mrflinray))) and (77)  (S ≦ p_(minslowmu))) {   v_(mrfrefy) = v_(mrflinray); }

A further magnitude limitation is necessary

if (v_(mrfrefy) > 0)   v_(mrfrefy) = min(min(p_(betabound)|v_(mrfx) |, p_(maxlatvel))),v_(mrfrefy)); else (78)   v_(mrfrefy) = max(max(−p_(betabound)|v_(mrfx) |,−p_(maxlatvel))),v_(mrfrefy));

If the conditions in (75) are not satisfied, the reference lateral velocity will be switched to the linear lateral velocity

else   {     v_(mrfrefy) = v_(mrflinray); (79)   }

Now the low frequency portion v_(mrfssy) of the lateral velocity at the rear axle and on the moving road frame in step 714 may be determined. This may be done through a steady state recovery filter in the following

v _(mrfssy) =d _(ssc) v _(mrfssy) +n _(ssc) v _(mrfrefy)   (80)

where d_(ssc),n_(ssc) are the filter coefficients.

Another method to obtain the low frequency portion v_(mrfssy) is through the following initial condition resetting

 v_(mrfssy) = d_(ssc)v_(mrfssy) + n_(ssc)v_(mrfrefy) if (( v_(mrfray)v_(mrflinry) < 0   or (v_(mrflinry) ≧ 0 and v_(mrfray) ≦ v_(mrflinry))   or (v_(mrflinry) ≦ 0 and v_(mrfray) ≧ v_(mrflinry)) (81)   ) and (S ≦ p_(minslowmu))) {   v_(mrfssy) = v_(mrflinry) − v_(mrfdyny); }

Now the computation of the dynamic portion of the lateral velocity at the rear axle on the moving road frame in step 716 is considered. This may be referred to as the high frequency portion. The derivative of such a velocity may be expressed as in the following

dv _(mrfy) =ā _(mrfy)−ω_(mrfz) v _(mrfx) v _(mrfz)   (82)

Considering the vertical velocity is small enough, the following is true

dv _(mrfy) =ā _(mrfy)−ω_(mrfz) v _(mrfx)   (83)

Using a anti-drift integration, the dynamic portion v_(mrfdyny) of the lateral velocity v_(mrfray) may be calculated as in the following

v _(mrfdyny) =d _(adi) v _(mrfdyny) +n _(adi) dv _(mrfy)   (84)

Since the potential errors in the computation in (83), a computation combining (84) with initial condition enforcement is conducted as in the following

v_(mrfdyny) = d_(adi)v_(mrfdyny) + n_(adi)dv_(mrfray) if (( v_(mrf)v_(mrflinry) < 0    or (v_(mrflinry) ≧ 0 and v_(mrfray) ≦ v_(mrflinry))    or (v_(mrflinry) ≦ 0 and v_(mrfray) ≧ v_(mrflinry)) (85)   ) and (S ≦ p_(minslowmu))) {   v_(mrfdyny) = v_(mrflinry) − v_(mrfssy); }

In step 718 the final lateral velocity may be determined by combining the steady state and dynamic lateral velocity. A variation of (85) may be further conducted by using a smooth ramping to replace a sudden switch

v_(mrfdyny) = d_(adi)v_(mrfdyny) + n_(adi)dv_(mrfray) if (( v_(mrf)v_(mrflinry) < 0    or (v_(mrflinry) ≧ 0 and v_(mrf) ≦ v_(mrflinry))    or (v_(mrflinry) ≦ 0 and v_(mrf) ≧ v_(mrflinry))   ) and (S ≦ p_(minslowmu))) (86) {   v_(mrfdyny) = v_(mrflinry) − v_(mrfssy) + (v_(mrfdyny) − v_(mrflinry) + v_(mrfssy))/ p_(rho); }

A further variation is the following

if (( v_(mrf)v_(mrflinry) < 0    or (v_(mrflinry) ≧ 0 and v_(mrf) ≦ v_(mrflinry))    or (v_(mrflinry) ≦ 0 and v_(mrf) ≧ v_(mrflinry))   ) and (S ≦ p_(minslowmu))) (87)   v_(mrfdyny) = v_(mrflinry) − v_(mrfssy) + (v_(mrfdyny) − v_(mrflinry) + v_(mrfssy))/ p_(rho); else   v_(mrfdyny) = d_(adi)v_(mrfdyny) + n_(adi)dv_(mrfray);

The final lateral velocity at the rear axle and on the moving road frame may be used to control various vehicle systems including but not limited to a yaw control system and other dynamic control systems.

Longitudinal Velocity Determination

In this embodiment, the vehicle longitudinal velocity is determined based on the compensated wheel speed sensor signals, the IMU sensor cluster outputs and the computed signals from the other modules in ISS unit, the wheel control status (in ABS, TCS, RSC, AYSC, or driver braking), engine torque and throttle information. The vehicle longitudinal velocity considered here is defined along the longitudinal direction of the moving road frame. The wheel speeds are the linear velocities which are projected on the instantaneous plane defined by the moving road frame (the longitudinal and lateral axle of the moving road frame). The algorithms and methods presented here reduce errors in existing longitudinal velocity computations used in brake control systems (ABS and TCS) and the vehicle stability control systems (such as ESC and RSC), which are based on the traditional ESC sensor set or the RSC sensor set.

The vehicle longitudinal velocity is an important variable used for vehicle dynamics controls. For example, the ABS and TSC brake controls activate the individual brake based on the longitudinal slip ratio of the corresponding wheel. The longitudinal slip ratio at a specific corner measures the relative difference between the velocity of the vehicle body at the corresponding corner and the linear velocity of the center of the corresponding wheel. Slip ratios are also used in lSCS functions such as RSC, YSC and LSC.

Another usage of the longitudinal velocity is to calculate the vehicle sideslip angle which is obtained by dividing the calculated lateral velocity at specific location of the vehicle body (for example at the vehicle center of gravity or at the middle of the rear axle) by the longitudinal velocity The computation of lateral velocity is affected by the accuracy of the longitudinal velocity, the errors in the longitudinal velocity itself will be directly transferred to the errors in sideslip angle computations

There exist many methods pursued by the existing stability controls, none of them achieves the accuracy and robustness that is deemed necessary for the ISCS function requirement. For example, in a driving scenario for a RWD vehicle equipped with the RSC system, during a two-wheel-lift event (two inside wheels are up in the air), the RSC requests a large brake pressure in the front outside wheel and there is some driving torque applied to the rear wheel. In this case, none of the four wheel speed sensor signals provide useful information to specify the vehicle longitudinal velocity. A natural approach is to integrate the available longitudinal acceleration to obtain an indication of the vehicle speed. However, due to the gravity contamination through vehicle pitch attitude and the lateral motion contamination due to the vehicle's lateral sliding motion, such an integration of the longitudinal acceleration may have significant errors and may further reduce the effectiveness of the RSC function or introduce potential false activations.

Another problem with the existing method is that on ice or snow the four individual wheels may all be in abnormal status (none of them is a good indication of the vehicle speed due to the fact that none of them is in pure rolling), and the longitudinal acceleration is again contaminated by the side sliding of the vehicle, the existing method generates significant errors which will affect the control actions.

Referring now to FIG. 9, the longitudinal velocity estimation is determined using a wheel speed block 800 that receives the wheel speeds and provides a scaling factor 802 and wheel compensation 804 therein. The output of the wheels' speed compensations is provided to a longitudinal reference, contact patch velocity monitoring and the longitudinal slip ratio block 806. The output of the longitudinal velocity is fed back into the scaling factor block 802 from block 806. Block 806 receives various inputs from sensors and the like. For example, the moving reference frame yaw rate, the steering angle, the wheels and powertrain actuation status, the driver braking status, and the linear sideslip angle at the moving road frame is provided to block 806. A longitudinal velocity derivative is determined using the longitudinal acceleration at the rear axle in the moving road frame, the yaw rate in the moving road frame, the G pitch and R pitch factors, and the lateral velocity at the rear angle in the moving road frame. The longitudinal velocity derivative is determined in block 808 from the previous four sensors and is provided to block 806 to determine a longitudinal slip ratio in the contact patch velocity monitor. The output of the block 806 is provided to arbitration and filtering block 310 which ultimately determines the longitudinal velocity.

Referring now to FIG. 10, the block 800 is illustrated in further detail. Block 800 includes the front left wheel speed-based algorithm 820, the front right wheel-based algorithm 822, the rear left wheel-based algorithm 824, and the rear right wheel-based algorithm 826. A front wheel-based algorithm 828 and a rear wheel-based algorithm 830 are also provided. Contact patch velocity estimation 832 and a slip ratio computation 834 may also be provided. Computations with these signals and the wheel actuation status and driver braking status may be provided in block 806. Block 840 determines whether the wheels are divergent. In block 842 if the wheels are divergent, the wheel speeds may be trusted and the longitudinal reference velocity is output therefrom. If the wheels are not divergent in block 842, an SSC filter may be provided. SSC filter combines the longitudinal velocity derivative with an anti-integration drift filter 846 to provide an output therefrom. The specifics of the above-mentioned calculations are further described below with respect to block 14.

Referring now to FIG. 11, in step 900 the various wheel speed sensors are determined. In step 902, the other sensor outputs are determined and translated into moving road plane.

The rate change of the vehicle longitudinal velocity in the moving road frame may be calculated from the signals in the moving road frame in step 904 in the following:

{dot over (v)} _(mrfx) =a _(mrfx)+ω_(mrfz) v _(mrfray) +g sin θ_(mrfy)   (88)

where θ_(mrfy) is the pitch angle of the moving road frame with respect to sea level, it is equivalent to the so-called road inclination angle or road slope. θ_(mrfy) may be calculated in the following

θ_(mrfy)=θ_(sy)−Δθ_(y)−θ_(yr)   (89)

where θ_(ry) is the global pitch angle of the sensor frame calculated based on the IMU sensor signals, Δθ_(y) is the sensor pitch misalignment and θ_(yr) is the relative pitch that reflects the suspension height difference between the front and the rear suspensions.

Since computing the longitudinal velocity at the location of the center of the rear axle is desired, the above computation in (88) will be further changed to

dv _(mrfrax) =a _(mrfr) +p _(xsl2ra)(ω_(mrfy) ²+ω_(mrfz) ²)+ω_(mrfz) v _(mrfray) +g sin θ_(mrfy)   (90)

where p_(xsl2ra) is the longitudinal distance between the sensor location and the rear axle and v_(mrfray) is the lateral velocity at the rear axle.

Theoretically, the longitudinal velocity at rear axle v_(mrfrax) is ready to be obtained from (90) through integration. However, due to sensor errors reflected in a_(mrfx)) ω_(mrfy) and ω_(mrfz), and the computational errors reflected in v_(mrfray) and θ_(mrfy), such a direct long-term integration will generate significant errors. Instead of relying heavily on the integration of dv_(mrfrax) all the time, the wheel speed sensor signals are used. A method to compute the longitudinal velocity which combines dv_(mrfrax) and wheel speeds is determined below.

In step 906, the wheel speeds may be used to determine the longitudinal velocity. When the wheel speed sensor signals are properly compensated with the correct scaling factors, and when the vehicle is under normal driving conditions, the compensated wheel speed sensor signals are the correct indications of the vehicle longitudinal velocity. In such driving cases, the wheels are usually operating in a stable condition, i.e., their motion is close to the free rolling. In other driving situations, when the wheels are far away from free rolling, wheel speeds are no longer the good indications of the vehicle longitudinal speeds. Such cases could happen when the wheels are braked due to a driver's braking request: the wheels are braked due to brake control requests, the wheels are driven with large driving torques (for driving wheels, which usually causes wheels' spinning), the wheels lose traction (driving torque applied to wheels are greater than the road traction limit), the wheels slide longitudinally, the wheels slide laterally (vehicle builds up large sideslip angle), or the wheels are up in the air (for example, the inside wheels during a sharp turn).

However if there exists at least one wheel in free rolling motion and the rest of wheels are in unstable conditions, the wheel speed may still be used for determining the vehicle longitudinal velocity.

If the front left wheel is in free rolling or close to free rolling, the vehicle longitudinal velocity at the middle point of the rear axle but projected on the moving road frame may be calculated based on this wheel's compensated wheel speed signal as in the following

$\begin{matrix} {{v(0)} = \frac{{\kappa_{0}w_{0}} + {\omega_{mrfz}\left\lbrack {{t_{f}{\cos (\delta)}} - {b\; {\sin (\delta)}}} \right\rbrack}}{{\cos (\delta)} + {{\sin (\beta)}{\sin (\delta)}}}} & (91) \end{matrix}$

where δ is the wheel's steering angle (road wheel steering angle), β is the sideslip angle at the rear axle, b is the vehicle base, t_(f) is the half track of the front axle; ω_(mrfz) is the vehicle body's yaw rate projected in the moving road frame, κ_(O) is the wheel speed scaling factor of the front left wheel.

Similarly, if the front right wheel is in free rolling or close to free rolling, the vehicle longitudinal velocity at the middle point of the rear axle but projected in the moving road frame may be calculated based on this wheel's compensated wheel speed signals as in the following

$\begin{matrix} {{v(1)} = \frac{{\kappa_{1}w_{1}} + {\omega_{mrfz}\left\lbrack {{t_{f}{\cos (\delta)}} - {b\; {\sin (\delta)}}} \right\rbrack}}{{\cos (\delta)} + {{\sin (\beta)}{\sin (\delta)}}}} & (92) \end{matrix}$

where κ₁ is the wheel speed scaling factor of the front right wheel.

If the rear left wheel is in free rolling, the vehicle longitudinal velocity at the middle point of the rear axle but projected in the moving road frame may be calculated based on this wheel's compensated wheel speed signals as in the following

v(2)=κ₂ω₂+ω_(mrfz) t _(r)   (93)

where t_(r) is the half track of the rear axle, κ₂ is the wheel speed scaling factor of the rear left wheel. If the rear left wheel is in free rolling, the vehicle longitudinal velocity at the middle point of the rear axle but projected in the moving road frame may be calculated based on this wheel's compensated wheel speed signals as in the following

v(3)=κ₃ω₃−ω_(mrfz) t _(r)   (94)

where κ₃ is the wheel speed scaling factor of the rear right wheel.

Based on the above four computed longitudinal velocities from the four wheel speed sensor signals, the velocity variables with their relative positions among the four variables may be computed as in the following

v _(f min)=min(v(0),v(1);v _(fav)=(v(0)+v(1))/2;v _(f max)=max(v(0),v(1));

v _(r min)=min(v(2),v(3));v _(rav)=(v(2)+v(3))/2;v _(r max)=max(v(2),v(3));

V _(min)=min(v _(f min) ,v _(r min));v _(max)=max(v _(f max) ,v _(rmax));v _(mid) =v _(fav) +v _(rav)−(v _(min) +v _(max))/2;   (95)

For a given vehicle longitudinal velocity v_(mrfrax), the linear velocity at each of the wheels but projected along the wheel longitudinal directions in step 908 can be computed as

v _(vmrfrax)(0)=v _(mrfrax) cos(δ)+ω_(mrfz)(b sin(δ)+t _(f) cos(δ))

v _(vmrfrax)(1)=v _(mrfrax) cos(δ)+ω_(mrfz)(b sin(δ)−t _(f) cos(δ))

v _(vmrfrax)(2)=v _(mrfrax)+ω_(mrfz) t _(r)

v _(vmrfrax)(3)=v _(mrfrax)−ω_(mrfz) t _(r)   (96)

Based on those linear corner speeds, the longitudinal slip ratios and their related quantities can be calculated in step 910 as in the following, where it is evident ρ(i) is the traditional definition of the longitudinal slip ratio of the ith wheel

for(i = 0; i < 4; i++) { ${{\rho (i)} = {\frac{v(i)}{\max \; \left( {v_{\min},{v_{vmrfrax}(i)}} \right)} - 1}};$ ρ(i) = min(1, max(−1, ρ(i)); ${{d\; {\rho (i)}} = {{p_{dlpf}d\; {\rho (i)}} + {\left( {1 - p_{dlpf}} \right)\frac{{\rho (i)} - {\rho_{z\; 1}(i)}}{\Delta T}}}};$ (97) ρ_(z1)(i) = ρ(i); ρ_(p)(i) = ρ(i)d ρ(i); ρ_(r)(i) = v_(vmrfrax) (i)d ρ(i);  

Now the quality of the wheel speed signals on which the vehicle longitudinal velocity can be calculated in step 912 is determined for each wheel as in the following

Also notice that the braking action at each of the wheels will affect using that wheel's wheel speed sensor to compute the longitudinal velocity. For this reason, the ith wheel's slip trend (divergent or convergent longitudinal slip indicated through the Boolean variable S[i].convergent_slip), acceleration trend (decelerated or accelerated wheel, indicated through the Boolean variable S[i].whldec and S[i].whlacc) and braking action (indicated through the Boolean variable S[i].inbraking) may be determined in the following

for(i = 0; i < 4; i++) { if (ρ_(p)(i) ≦ 0 ) S[i].convergent_slip = 1;  else S[i].convergent_slip = 0;  if ( ρ_(r)(i) > p_(minwacc)) S[i].whlacc = 1;  else if (ρ_(r)(i) < −p_(minwacc)) S[i].whldec = 1;  else (98)  { S[i].whlacc = 0;  S[i].whldec = 0; }   if (S[i].ABS = 0 & &S[i].BTCS = 0 & & S[i].YSC = 0     & & S[i].RSC =0 &&S[i].LSC = 0 & & P(i) ≦ p_(minpres))    S[i].inbraking = 0;   else   S[i].inbraking = 1;   } where S[i].ABS , S[i].BTCS , S[i].YSC , S[i].RSC and S[i].LSC are Boolean variables whose active values mean there are ABS, BTCS (brake traction control), YSC, RSC and LSC brake activations respectively.

The contact patch velocity difference between the left and right wheels at the same axles and the average of the contact patch velocities at the same axles can be computed as in the following

$\begin{matrix} {{{{{\Delta \; v_{cpfa}} = {{\frac{w_{1} - w_{0}}{2} - {\omega_{mrfz}t_{f}{\cos (\delta)}}}}};}{{v_{cpfa} = {{\frac{w_{1} + w_{0}}{2} - {\omega_{mrfz}b\; \sin \; (\delta)} - {v_{mrfrax}\left( {{\cos (\delta)} + {{\sin (\delta)}\beta_{mrfra}}} \right)}}}};}{{{\Delta \; v_{cpra}} = {{\frac{w_{3} - w_{2}}{2} - {\omega_{mrfz}t_{r}}}}};}\; {v_{cpra} = {{\frac{w_{3} + w_{2}}{2} - v_{mrfrax}}}}};} & (99) \end{matrix}$

and they can be used to determine if the contact patches are drastically different from each other at the same axle. Such a logic is conducted as in the following

if (Δv_(cpfa) ≦ p_(dfth)v_(mrffavx)) & &v_(cpfa) ≦ p_(mfth)v_(mrffavx)) { S[0].convergent_cpv = 1;S[1].convergent_cpv = 1; } else { S[0].convergent_cpv = 0;S[1].convergent_cpv = 0; } if (Δv_(cpra) ≦ p_(dfth)v_(mrfravx) & &v_(cpra) ≦ p_(mfth)v_(mrfravx)) (100) { S[2].convergent_cpv = 1;S[3].convergent_cpv = 1; } else { S[2].convergent_cpv = 0;S[3].convergent_cpv = 0; }

Next the drive train mode and torques are determined in step 914 based on the driving torque distribution D_(r) information

if (D_(τ) ≦ 0)   tDriveMode = REAR; else if (D_(τ) > 0 & & D_(τ) < p_(lrgdistribution))   tDriveMode = TOD; (101) else   tDriveMode = FOUR; where REAR means RWD, TOD means torque-on-demand, and FOUR means ×4 mode. p_(irgdistribution) is a threshold for torque distribution to determine the drive mode. Notice that a more detailed drivetrain determination can be used.

The braking torque applied to each of the four wheels can be computed based on the caliper pressure P(i) at the ith wheel

τ_(b)(0)=P(0)*p _(brkTGain0);τ_(b)(1)=P(1)*p _(brkTGain1);

τ_(b)(2)=P(2)*p _(brkTGain2);τ_(b)(3)=P(3)*p _(brkTGain3);   (102)

where p_(brkTGain) is the braking torque gain for the ith wheel.

Reference Velocity by Integrating Wheel Speed Based Computation With the Acceleration Based Computation

Based on the driver mode and the computed braking torque, the total torque applied to the ith wheel can be calculated. For example, in the following, for RWD and FWD modes, the total wheel end torque can be calculated using the engine torque distributed to the axle, τ_(axletorque)

if (tDriveMode = REAR) { τ_(d)(0) = (−1)*τ_(b)(0); τ_(d)(1) = (−1)*τ_(b)(1); $\left. {{{\tau_{d}(2)} = {\frac{\tau_{axletorque}}{2} - {\tau_{b}(2)}}};{{\tau_{d}(3)} = {\frac{\tau_{axletorque}}{2} - {\tau_{b}(3)}}};}\mspace{25mu} \right\}$ (103) if (tDriveMode = FRONT) { ${{{\tau_{d}(0)} = {\frac{\tau_{axletorque}}{2} - {\tau_{b}(0)}}};{{\tau_{d}(1)} = {\frac{\tau_{axletorque}}{2} - {\tau_{b}(1)}}};}\;$ τd(2) = (−1)*τ_(b)(2); τ_(d)(3) = (−1)*τ_(b)(3);  }

In step 916, the reference longitudinal velocity is used to find the low frequency portion of the vehicle longitudinal velocities in the following.

Non-wheel lifting case which is indicated by relative roll, wheel departure angle and the front wheel control bits

f (| θ_(xr) |≦ p_(75th)p_(rollgradient) & & (θ_(wda) = 0) & & S[0].RSC = 0 & &S[1].RSC = 0 ) {  if (tDriveMode = REAR) (104)  {

Considering the torque-on-demand drivetrain TOD is a transitional mode, a 4×2 assumption for TOD may be used. For some other powertrains, TOD may not be a transitional mode.

For the four-wheel-coasting-case, the reference longitudinal velocity is ramping to the maximum longitudinal velocity v_(max)

f (  S[0].inbraking = 0 & &S[1].inbraking = 0 & &S[2].inbraking = 0 & &S[3].inbraking = 0 & &Θ_(throttle) ≦ 0 & &τ_(axletorque) ≦ p_(maxcoasttorq) ) (105) { $\left. {{v_{mrfraxref} = {v_{\max} + \frac{\left( {v_{mrfraxref} - v_{\max}} \right)}{\xi}}};}\mspace{31mu} \right\}$ where Θ_(throttle) is the throttle angle, J is a constant greater than 1 which is related to the ramping rate, p_(maxcoasttorq) is a threshold which is the maximum allowed torque for the vehicle to be coasted.

If the above four-wheel-coasting-conditions are not satisfied, the following two-front-wheel-coasting case is performed (one or two rear wheels in braking, or rear wheels in driving):

else if (S[0].inbraking = 0 & &S[1].inbraking = 0) { if (S[0].within_smallband = 1) ${v_{mrfraxref} = {{v(0)} + \frac{\left( {v_{mrfraxref} - {v(0)}} \right)}{\xi}}};$ else if (S[1].within_smallband = 1) ${v_{mrfraxref} = {{v(1)} + \frac{\left( {v_{mrfraxref} - {v(1)}} \right)}{\xi}}};$ else if (Θ_(throttle) ≦ 0 & &τ_(axletorque) ≦ p_(maxtorquecost)) { if (v_(fmax) > p_(minvel)) ${v_{mrfraxref} = {v_{f\; \max} + \frac{\left( {v_{mrfraxref} - v_{f\; \max}} \right)}{\xi}}};$ else  v_(mrfraxref) += dv_(mrfrax)ΔT;  } else if ( Θ_(throttle) > 0 ∥ τ_(axletorque) > p_(maxtorquecost)) $\left. \left. {{v_{mrfraxref} = {v_{fo} + \frac{\left( {v_{mrfraxref} - v_{fo}} \right)}{\xi}}};}\mspace{31mu} \right\} \mspace{40mu} \right\}$ (106)

If the four-wheel-coasting-condition and the two-wheel-coasting conditions are not satisfied, the two-rear-wheel coasting case is checked (one or two front wheels in braking):

else if (S[2].inbraking = 0 & &S[3].inbraking = 0 & &Θ_(throttle) > 0 && τ_(axletorque) > p_(maxtorquecost)) { if (S[2].within_smallband = 1) ${v_{mrfraxref} = {{v(2)} + \frac{\left( {v_{mrfraxref} - {v(2)}} \right)}{\xi}}};$ else if (S[3].within_smallband = 1) ${v_{mrfraxref} = {{v(3)} + \frac{\left( {v_{mrfraxref} - {v(3)}} \right)}{\xi}}};$ else if (S[0].inbraking = 1 & &S[1].inbraking = 1) { if (v_(fmax) > p_(minvel) & & |ρ(i_(m)) |≦ p_(slipratiobd)) ${v_{mrfraxref} = {v_{r\; \max} + \frac{\left( {v_{mrfraxref} - v_{r\; \max}} \right)}{\xi}}};$ else  v_(mrfraxref) += dv_(mrfrax)ΔT;  } else $\left. {{v_{mrfraxref} = {v_{rav} + \frac{\left( {v_{mrfraxref} - v_{rav}} \right)}{\xi}}};}\mspace{31mu} \right\}$ (107)

If the above four-wheel-coasting and two-wheel-coasting conditions are not satisfied, the single wheel coasting case for the front-left wheel is checked:

else if (S[0].within_middle = 1 & &S[0].convergent_slip = 1 & &S[0].convergent_cpv = 1 & &S[0].inbraking = 0 & &N(0) > p_(fstatical)(1 − p_(nlloss))) (108) { $\left. {{v_{mrfraxref} = {{v(0)} + \frac{\left( {v_{mrfraxref} - {v(0)}} \right)}{\xi}}};}\mspace{31mu} \right\}$

If the four-wheel-coasting, two-wheel-coasting and front-left single-wheel-coasting conditions are all not satisfied, the front right wheel coasting condition is checked:

else if (S[1].within_middle = 1 & &S[1].convergent_slip = 1 & &S[1].convergent_cpv = 1 & &S[1].inbraking = 0 & &N(1) > p_(fstatical)(1 − p_(nlloss))) (109) { $\left. {{v_{mrfraxref} = {{v(1)} + \frac{\left( {v_{mrfraxref} - {v(1)}} \right)}{\xi}}};}\mspace{31mu} \right\}$

If all the above checked conditions are not satisfied, the rear-right wheel for coasting conditions is then checked:

else if (S[2].within_middle = 1 & &S[2].convergent_slip = 1 & &S[2].convergent_cpv = 1 & &S[2].inbraking = 0 & &N(2) > p_(fstatical)(1 − p_(nlloss))&&Θ_(throttle) ≦ 0 (110) &&τ_(axletorque) > p_(maxtorquecoast)&&|β_(mrfra) |≦ p_(ssamidbd)) { $\left. {{v_{mrfraxref} = {{v(2)} + \frac{\left( {v_{mrfraxref} - {v(2)}} \right)}{\xi}}};}\mspace{31mu} \right\}$

If none of the four-wheel-coasting, two-wheel-coasting, and front-left, front-right and rear-left single wheel coasting conditions are satisfied, the front right wheel for coasting is then checked:

else if (S[3].within_middle = 1 & &S[3].convergent_slip = 1 & &S[3].convergent_cpv = 1 & &S[3].inbraking = 0 & &N(3) > p_(fstatical)(1 − p_(nlloss))&&Θ_(throttle) ≦ 0 (111) &&τ_(axletorque) > p_(maxtorquecoast)&&|β_(mrfra) |≦ p_(ssamidbd)) { $\left. {{v_{mrfraxref} = {{v(3)} + \frac{\left( {v_{mrfraxref} - {v(3)}} \right)}{\xi}}};}\mspace{31mu} \right\}$ }

If the vehicle is not in 4×2 mode (BWD mode as discussed before, FWD can be similarly considered), but in 4×4 mode (that is, in all-wheel-drive or toque-on-demand or 4-high or 4-low modes), the following check is performed for a single wheel condition regardless what the other three wheels are doing (could be braking, driving or coasting)

else { if (S[0].within_smallband = 1) ${v_{mrfraxref} = {{v(0)} + \frac{v_{mrfraxref} - {v(0)}}{\xi}}};$ else if (S[1].within_smallband = 1) ${v_{mrfraxref} = {{v(1)} + \frac{v_{mrfraxref} - {v(1)}}{\xi}}};$ (112) else if ( S[2].within_smallband = 1) ${v_{mrfraxref} = {{v(2)} + \frac{v_{mrfraxref} - {v(2)}}{\xi}}};$ else if (S[3].within_smallband = 1) ${v_{mrfraxref} = {{v(3)} + \frac{v_{mrfraxref} - {v(3)}}{\xi}}};$

If none of the above single wheel condition is satisfied, the following convergent wheel speed case and the four wheel coasting case is used:

else if (|v_(max) − v_(min) |< p_(smith) | v_(max) + v_(min) |) ${v_{mrfraxref} = {v_{mid} + \frac{v_{mrfraxref} - v_{mid}}{\xi}}};$ else if (S[0].inbraking = 0 & &S[1].inbraking = 0 & &S[2].inbraking = 0 & &S[3].inbraking = 0 (113) &&(Θ_(throttle) ≦ 0 & &τ_(axletorque) ≦ p_(maxtorquecoast))) { $\left. {{v_{mrfraxref} = {{\left( {v_{\max} + v_{\min}} \right)\text{/}2} + \frac{v_{mrfraxref} - {\left( {v_{\max} + v_{\min}} \right)\text{/}2}}{\xi}}};}\mspace{31mu} \right\}$

When the vehicle is driven with throttle off, the following computations are conducted:

else if (Θ_(throttle) ≦ 0 & &τ_(axletorque) ≦ p_(maxtorquecoast))) { if (S[0].within_smallband = 1) ${v_{mrfraxref} = {{v(0)} + \frac{v_{mrfraxref} - {v(0)}}{\xi}}};$ else if ( S[1].within_smallband = 1) $v_{mrfraxref} = {{v(1)} + \frac{v_{mrfraxref} - {v(1)}}{\xi}}$ else if (S[2].within_smallband = 1) $v_{mrfraxref} = {{v(2)} + \frac{v_{mrfraxref} - {v(2)}}{\xi}}$ else if (S[3].within_smallband = 1) ${v_{mrfraxref} = {{v(3)} + \frac{v_{mrfraxref} - {v(3)}}{\xi}}};$ else if (| v_(max) − v_(min) |< p_(smith) | v_(max) − v_(min) |) ) $\left. \mspace{31mu} {{{v_{mrfraxref} = {v_{mid} + \frac{v_{mrfraxref} - v_{mid}}{\xi}}};}{{{{else}\mspace{14mu} v_{mrfraxref}} = {v_{\max} + \frac{v_{mrfraxref} - v_{\max}}{\xi}}};}}\mspace{31mu} \right\}$ (114)

If there are no braking on any of the wheels, then the following is performed:

else if (S[0].inbraking = 0 & &S[1].inbraking = 0 & &S[2].inbraking = 0 & &S[3].inbraking = 0 ) { if ( S[0].within_smallband = 1) ${v_{mrfraxref} = {{v(0)} + \frac{v_{mrfraxref} - {v(0)}}{\xi}}};$ else if ( S[1].within_smallband = 1) ${v_{mrfraxref} = {{v(1)} + \frac{v_{mrfraxref} - {v(1)}}{\xi}}};$ else if (S[2].within_smallband = 1) ${v_{mrfraxref} = {{v(2)} + \frac{v_{mrfraxref} - {v(2)}}{\xi}}};$ else if (S[3].within_smallband = 1) ${v_{mrfraxref} = {{v(3)} + \frac{v_{mrfraxref} - {v(3)}}{\xi}}};$ else if (| v_(max) − v_(min) |<| v_(max) + v_(min) | p_(smith)) ${v_{mrfraxref} = {v_{mid} + \frac{v_{mrfraxref} - v_{mid}}{\xi}}};$ else   $\left. {{v_{mrfraxref} = {v_{\min} + \frac{v_{mrfraxref} - v_{\min}}{\xi}}};}\mspace{31mu} \right\}$ (115)

If none of the above conditions is satisfied, the integration of the time derivative of the longitudinal velocity along the moving road frame is conducted:

    else       v_(mrfraxref) + = dv_(mrfrax) * ΔT;      } (116)   } }

For a potential inside wheel/wheel lifting case:

else { if ( Θ_(throttle) ≦ 0 & &τ_(axletorque) ≦ p_(maxtorquecoast) & &dv_(mrfraxf) ≦ 0 & &dv_(mrfraxf) >−p_(mindeepbrakedvc) & &|ρ(i_(ro)) |≦ p_(slipratiobd) { $\left. {{v_{mrfraxref} = {v_{ro} + \frac{v_{mrfraxref} - v_{ro}}{\xi}}};}\mspace{31mu} \right\}$ else if ( (τ_(axletorque) ≧ p_(mintorquehmu) ∥ Θ_(throttle) ≧ p_(minthrottlehmu)) (117) & &dv_(mrfraxf) > 0 & &dv_(mrfraxf) ≦ p_(mindeepbrakedvc)) ${v_{mrfraxref} = {v_{fmax} + \frac{v_{mrfraxref} - v_{fmax}}{\xi}}};$ else  v_(mrfraxref) += dv_(mrfrax)*ΔT; } } }

In step 918, The longitudinal velocity may then be computed based on the reference longitudinal velocity and a washout integration of the derivative of the longitudinal velocity such that the non-smooth switch in the reference longitudinal velocity is

if (v_(xmrf) ≧ p_(smspd)) { if ((v_(max) − v_(min)) > p_(smpe) (v_(max) + v_(min))) { //Compute High Frequency Portion of Longitudinal Velocity v_(mrfxdyn) = p_(adid)v_(mrfxdyn) + p_(adin)dv_(mrfrax); //Compute Low Frequency Portion of Longitudinal Velocity v_(mrfxss) = p_(adid)v_(mrfxss) + p_(sscn)v_(mrfraxref); v_(mrfx) = v_(mrfxdyn) + v_(mrfxss); v_(mrfxls) = v_(mrfx);   v_(mrfyls) = v_(mrfy);  } //Compute Longitudinal Velocity During Convergent Wheel Speeds else if ((v_(max) − v_(min)) ≦ p_(smpc) (v_(max) + v_(min))) { ${v_{mrfx} = {v_{mid} + \frac{v_{mrfx} - v_{mid}}{\xi}}};$ ${V_{mrfxss} = {V_{mrfx} + \frac{V_{mrfxss} - V_{mrfx}}{\xi}}};\; {V_{mrfxdyn} = {V_{mrfx} - V_{mrfxss}}};$ v_(mrfxls) = v_(mrfx);  v_(mrfyls) = v_(mrfy);  } (118) }

If during low vehicle speed,

else { //Compute Longitudinal Velocity During Low Speed Driving v_(mrfyls) + = dv_(mrfray)ΔT; dv_(mrfxls) = a_(mrfxaug) + w_(mrfz)v_(mrfyls); v_(mrfxls) + = dv_(mrfrax)ΔT; if (v_(mrfxls) < 0 & &v_(mrfxls) ≧ −v_(max))    v_(mrfxls) = −v_(max); (119)  if (v_(mrfxls) ≧ p_(smspd))   v_(mrfxls) = p_(smspd);  v_(mrfxdyn) = p_(adid)v_(mrfxdyn) + p_(adin)dv_(mrfxls);  v_(mrfxss) = p_(adid)v_(mrfxss) + p_(sscn)v_(mrfxls);  v_(mrfx) = v_(mrfxdyn) + v_(mrfxss); }

The derivative of the longitudinal velocity based on the computed longitudinal velocity

dv _(mrfx) =p _(dervd1) dv _(mrfx−1) −p _(dervd2) dv _(mrfx 2) +p _(dervn)(v _(mrfx) −v _(mrfx−2));

dv _(mrfx−2) =dv _(mrfx−1);

dv _(mrfx−1) =dv _(mrfx);

v _(mrfx−2) =v _(mrfx−1);

v _(mrfx−1) =v _(mrfx);   (120)

The Wheel Speed Sensor Scaling Factor

In step 920, the wheel speeds scaling factors are determined such that the products of the scaling factors and the wheel speed sensor outputs provide the true linear speeds of the wheel centers during the wheels' coasting (free rolling). The final individual scaling factors involve several separated scaling factor computations, which include the dynamic scaling factors and several static scaling factors including: the left-to-right delta scaling factor at front axle; the left-to-right delta scaling factor at rear axle; the front-to-rear reference delta scaling factor; the global scaling factors.

The static scaling factors are calculated through several sequential computations,

In the first step, the left-to-right delta scaling factor at front axle or at the rear axle are first identified based on the individual longitudinal velocities calculated from the individual wheel speed sensor output, to try to balance the left and right wheel speed at the same axle during normal or costing wheel conditions.

In the second step, the reference wheel at front and the reference wheel at the rear are picked to compute the delta scaling factor, which is called the front-to-rear reference delta scaling factor.

In the third step, a reference wheel which does not need to have any scaling factor adjustment and which is one of the front and rear reference wheels used in step 2. The global scaling factor is calculated based on this reference wheel through matching its acceleration with the vehicle's longitudinal acceleration.

In the fourth step, the final static scaling factors are generated based on the delta scaling factors and the global scaling factors calculated in step 1 to 3.

In parallel to the static scaling factor determination, the dynamic scaling factors are determined based on the dynamic event.

The linear wheel speeds are used to describe the velocities of the centers of the wheels along their longitudinal directions. This information can be used to assess the longitudinal wheel slips used in brake and traction controls (such as ABS, TCS, RSC and ESC) and used to determine the longitudinal velocity of the vehicle. Such linear wheel speeds are translational velocities and cannot be directly measured unless expensive velocity sensors are used. The wheel speed sensor was invented for this purpose, however, it is rather misleading to call the output from a wheel speed sensor as the linear wheel speed since the wheel speed sensor really measures the wheels rotational speed and its output is the product of such rotational rate and the nominal rolling radius of the wheel. When the wheel's true rolling radius differs from the nominal value, the wheel speed sensor output will not be able to reflect the true vehicle velocity.

During the travel of a vehicle, many error sources contribute to the variation of the rolling radii and the following lists some of them: low inflation pressure in the tire of the interested wheel (the lower pressure, the smaller the rolling radius); uneven static loading distribution in the vehicle caused unevenly distributed rolling radii; the spare tire usage (the spare tire usually has a rolling radius which is about 15% to 17% smaller than the nominal rolling radius); the off-spec tire usage (customers sometimes use large off-spec tires, which have larger rolling radii); lateral load transfer during dynamic maneuvers such as potential RSC event (the load transfer may generate about 5% reduction in tire rolling radii for outside wheels during a one or two wheel lift event).

The above error sources will introduce errors in wheel linear speed determination, which further affects the accuracy of both wheel and the vehicle state estimations. For example, one of the important wheel states is the wheel's longitudinal slip ratio, its determination depends on the linear speed of the wheel and the longitudinal velocity of the vehicle. An erroneous computation in the wheel and vehicle states will ultimately cause false activations in controls or undermine the effectiveness of the control functions. Therefore it is desirable to counteract to the afore-mentioned error sources by adaptively adjusting the wheel scaling factors, which are the ratios between the true wheel rolling radii and the nominal wheel rolling radii, so as to adaptively calibrate the wheel speeds.

By using the nominal rolling radius, the wheel speed sensor outputs can be expressed as

w_(i)=ω_(is)R_(i)   (121)

for i=0,1,2,3 (indicating front-left, front-right, rear-left and rear-right wheels), where w_(i) is the well-known (but misleading) wheel speed sensor output for the ith wheel which depends on its rolling radius; cos is the rotational velocity of the ith wheel which is independent of its rolling radius; R_(i) is the nominal rolling radius of the ith wheel. Usually the front two wheels and the rear two wheels have the similar rolling radii

R₀=R₁=R_(f)

R₂=R₃=R_(r)   (122)

If the true rolling radius of the ith wheel is denoted as r_(i), then the ith wheel speed scaling factor is defined as

$\begin{matrix} {\chi_{i} = \frac{r_{i}}{R_{i}}} & (123) \end{matrix}$

Due to the computational advantage, instead of using (123), the delta scaling factors are used. The left-to-right delta scaling factors for the front wheels (Δχ_(f)) and rear wheels (Δχ_(r)) are defined as

$\begin{matrix} {{{\Delta \; \chi_{f}} = \frac{r_{0} - r_{1}}{R_{f}}}{{\Delta \; \chi_{r}} = \frac{r_{2} - r_{3}}{R_{r}}}} & (124) \end{matrix}$

The front-to-rear delta scaling factor is defined as

$\begin{matrix} {{\Delta \; \chi_{f\; 2r}} = \frac{2\left( {r_{ireff} - r_{irefr}} \right)}{R_{f} + R_{r}}} & (125) \end{matrix}$

where ireff is one of the front wheels which is called the front reference wheel, whose rolling radius does not need to be adjusted; irefr is one of the rear wheels which is called the rear reference wheel, whose rolling radius does not need to be adjusted.

The wheel scaling factor defined in (123) can also be determined though longitudinal acceleration alignment as in the following

$\begin{matrix} {\kappa_{g} = {{\frac{{\overset{.}{v}}_{x - {MRF}}}{{\overset{.}{v}}_{i}}\mspace{14mu} {if}\mspace{14mu} {\overset{.}{v}}_{i}} \neq 0}} & (126) \end{matrix}$

where {dot over (v)}_(x-MRF) is the time derivative of the longitudinal velocity obtained from the acceleration sensor signals and the calculated vehicle attitudes. Notice that {dot over (v)}_(x-MRF) is projected on the same plane on which the wheel travels. {dot over (v)}_(i) is the time derivative of the longitudinal velocity from the ith wheel speed sensor (but the turn effect, or say the difference due to vehicle turning is compensated).

Since the rolling radii are not known prior to the computation, (123)-(125) are not useful for the scaling factor computations. By using the sensed wheel angular velocities (wheel speed sensor outputs) and the longitudinal velocity at the center of the wheels, the rolling radii in (123)-(125) can be further expressed as

$\begin{matrix} {{{r_{0} = {\frac{v_{{wc}\; 0}}{w_{0}}R_{f}}};{r_{1} = {\frac{v_{{wc}\; 1}}{w_{1}}R_{f}}}}{{r_{2} = {\frac{v_{{wc}\; 2}}{w_{2}}R_{r}}};{r_{3} = {\frac{v_{{wc}\; 3}}{w_{3}}R_{r}}}}} & (127) \end{matrix}$

where v_(wci) is the linear longitudinal velocity of the center of the ith wheel. Therefore (124) can be expressed as

$\begin{matrix} {{{\Delta \; \chi_{f}} = {\frac{v_{{wc}\; 0}}{w_{0}} - \frac{v_{{wc}\; 1}}{w_{1}}}}{{\Delta \; \chi_{r}} = {\frac{v_{{wc}\; 2}}{w_{2}} - \frac{v_{{wc}\; 3}}{w_{3}}}}} & (128) \end{matrix}$

(128) is still not particularly useful due to the fact that v_(wci) are unknowns. Considering if the vehicle is driven on a straight road, ideally all v_(wci)s will be the same, hence if the following relative computations are introduced

$\begin{matrix} {{{\Delta\kappa}_{f} = \frac{r_{0} - r_{1}}{\left( {r_{0} + r_{1}} \right)/2}}{{\Delta\kappa}_{r} = \frac{r_{2} - r_{3}}{\left( {r_{2} + r_{3}} \right)/2}}{{\Delta\kappa}_{f\; 2r} = \frac{\left( {r_{i} - r_{j}} \right)}{\left( {r_{i} + r_{j}} \right)/2}}} & (129) \end{matrix}$

The unknowns v_(wci)s could be eliminated to obtain the following computations which only depend on the known variables

$\begin{matrix} {{{\Delta\kappa}_{f} = \frac{w_{1} - w_{0}}{\left( {w_{1} + w_{0}} \right)/2}}{{\Delta\kappa}_{r} = \frac{w_{3} - w_{2}}{\left( {w_{3} + w_{2}} \right)/2}}{{\Delta\kappa}_{f\; 2r} = \frac{\left( {w_{ireff} - w_{irefr}} \right)}{\left( {w_{ireff} + w_{irefr}} \right)/2}}} & (130) \end{matrix}$

Notice there are two issues associated with (130). The first issue is the contamination of the vehicle turning effect. That is, when a vehicle is driven in a turn, the outside wheels travel more distances in order to catch up with the inside wheels and they speed up, which would look much the same as the outside wheels have smaller rolling radii. However, if the corner wheel speed sensor outputs is transferred to the same location, this turning effect can be eliminated. The second issue is that (130) is clearly the determination of the relative scaling factors with respect to the average wheel speed of the right and the left wheels for Δκ_(f) and Δκ_(r), and with respect to the average wheel speed of the front and rear reference wheels for Δκ_(f2r).

The scaling factor difference of one wheel with respect to the other is important. A better definition for the front delta scaling factor would be

$\begin{matrix} {{\Delta\kappa}_{f} = \frac{w_{1} - w_{0}}{\max \left( {w_{1},w_{0}} \right)}} & (131) \end{matrix}$

However due to the later least computation method, if the noises in wheel speeds cause different wheels picked for the max, the above may not be robust. Instead of using (131) a correction to the computation (130) is conducted as in the following

$\begin{matrix} {{{{\Delta\kappa}_{f}^{\prime} = \frac{w_{1} - w_{0}}{\left( {w_{1} + w_{0}} \right)/2}},{{\Delta\kappa}_{f} = \frac{{\Delta\kappa}_{f}^{\prime}}{1 - {{{\Delta\kappa}_{f}^{\prime}}/2}}}}{{{\Delta\kappa}_{r}^{\prime} = \frac{w_{3} - w_{2}}{\left( {w_{3} + w_{2}} \right)/2}},{{\Delta\kappa}_{r} = \frac{{\Delta\kappa}_{r}^{\prime}}{1 - {{{\Delta\kappa}_{r}^{\prime}}/2}}}}{{{\Delta\kappa}_{f\; 2r}^{\prime} = \frac{\left( {w_{j} - w_{i}} \right)}{\left( {w_{j} + w_{i}} \right)/2}},{{\Delta\kappa}_{f\; 2r} = \frac{{\Delta\kappa}_{f\; 2r}^{\prime}}{1 - {{{\Delta\kappa}_{f\; 2r}^{\prime}}/2}}}}} & (132) \end{matrix}$

Overcoming the first issue by restricting the computation straight line driving not only limits the screening conditions for the computations but could also causes robustness issues, since a true straight line driving is hard to be determined. In order to solve the first issue, using the wheel speeds with removed turning effects is considered by compensating the yaw rate sensor signals. That is, the individual longitudinal velocities are worked with, which are calculated by transferring the wheel speed signals to the same location. In order to overcome the second issue, the acceleration calibration as in (136) will be used to compute portion of the scaling factor which could not be determined through Δκ_(f),Δκ_(r) and Δκ_(f2r).

An individual longitudinal velocity is defined as the one which measures the vehicle longitudinal velocity at the rear axle location which is calculated by transferring the four individual wheel speed sensor signals to same rear axle location. The following provides a computation of the individual longitudinal velocities

$\begin{matrix} {{v_{0} = \frac{w_{0} + {\omega_{mrfz}\left\lbrack {{t_{f}{\cos (\delta)}} - {b\; {\sin (\delta)}}} \right\rbrack}}{{\cos (\delta)} + {{\sin (\beta)}{\sin (\delta)}}}}{v_{1} = \frac{w_{1} - {\omega_{mrfz}\left\lbrack {{t_{f}{\cos (\delta)}} - {b\; {\sin (\delta)}}} \right\rbrack}}{{\cos (\delta)} + {{\sin (\beta)}{\sin (\delta)}}}}{v_{2} = {w_{2} + {ct}_{r}}}{v_{3} = {w_{3} - {\omega_{mrfz}t_{r}}}}} & (133) \end{matrix}$

where t_(f) and t_(r) are the half tracks for the front and rear axles, b is the vehicle base, ω_(mrfz) is the vehicle's yaw rate projected on the axle which is perpendicular to the plane on which the wheels are traveling, which is called moving road frame.

Based on the above individual longitudinal velocities, compute the delta scaling factors may be approximately computed as in the following

$\begin{matrix} {{{{\Delta\kappa}_{f}^{\prime} = \frac{v_{0} - v_{1}}{\left( {v_{0} + v_{1}} \right)/2}},{{\Delta\kappa}_{f} = \frac{{\Delta\kappa}_{f}^{\prime}}{1 - {{{\Delta\kappa}_{f}^{\prime}}/2}}}}{{{\Delta\kappa}_{r}^{\prime} = \frac{v_{2} - v_{3}}{\left( {v_{2} + v_{3}} \right)/2}},{{\Delta\kappa}_{r} = \frac{{\Delta\kappa}_{r}^{\prime}}{1 - {{{\Delta\kappa}_{r}^{\prime}}/2}}}}} & (134) \end{matrix}$

Using the computations from (134) to adjust the wheel speeds leads to the following

w _(i) ^(comp)=(1−|Δκ_(f)|)w _(i)

w _(j) ^(comp)=(1−|Δκ_(r)|)w _(j)   (135)

After adjustment in (135) is done, the compensated wheel speeds in step 922 will generate the same individual wheel longitudinal velocities in the front two wheels, or the same individual wheel longitudinal velocities in the rear two wheels.

Notice that the balanced individual longitudinal velocities could still be different from the true values due to the relative feature of the computation method in (134). Hence there is a need to do further compensation so as to balance the front and the rear individual longitudinal velocities.

Let ireff be a front reference wheel which is a wheel used as a reference to balance the front two wheels, and irefr be a rear reference wheel, which is used as a reference to balance the rear two wheels, then

$\begin{matrix} {{{\Delta\kappa}_{f\; 2r}^{\prime} = {2\frac{v_{ireff} - v_{irefr}}{v_{ireff} + v_{irefr}}}}{{\Delta\kappa}_{f\; 2r} = \frac{{\Delta\kappa}_{f\; 2r}^{\prime}}{1 - {{{\Delta\kappa}_{f\; 2r}^{\prime}}/2}}}} & (136) \end{matrix}$

Due to the noise and disturbance in the wheel speed sensor output, none of the computations in (134) and (136) have practical significance. In order to removing such noises and disturbances factors, a least square parameter identification method is used here. The least square parameter identification method is an averaging method which tries to find the optimal parameters to fit certain amount of data point. The least square algorithm for computing Δκ_(f) in (134) can be expressed as in the following

if (front wheel screening condition holds) { n_(f) = n_(f) + 1; if (n_(f) < p_(totlsstesps)) { A_(f) = A_(f) + (v₀ + v₁)²/4; B_(f) = B_(f) + (v₀ ² − v₁ ²)/2;  } else { ${{\Delta \; \kappa_{f}} = {\max \; \left( {{- p_{deviation}},{\min \; \left( {p_{deviation},\frac{B_{f}}{A_{f}}} \right)}} \right)}};$ ${{\Delta\kappa}_{f}^{\prime} = \frac{{\Delta\kappa}_{f}^{\prime}}{1 - {{{\Delta\kappa}_{f}^{\prime}}\text{/}2}}};$ A_(f) = 0; B_(f) = 0; N_(f) = 0; Δκ_(f) = p_(dlpf)Δκ_(f) + (1 − p_(dlpf))Δκ^(′) _(f);  } (137) }

The least square algorithm for computing Δκ_(r) in (134) can be expressed as in the following

if (rear wheel screening condition holds) { n_(r) = n_(r) + 1; if (n_(r) < p_(totlsstesps)) { A_(r) = A_(r) + (v₂ + v₃)²/4; B_(r) = B_(r) + (v_(r) ² − v₃ ²)/2;  } else { ${{\Delta \; \kappa_{r}} = {\max \; \left( {{- p_{deviation}},{\min \; \left( {p_{deviation},\frac{B_{r}}{A_{r}}} \right)}} \right)}};$ ${{\Delta\kappa}_{r}^{\prime} = \frac{{\Delta\kappa}_{r}^{\prime}}{1 - {{{\Delta\kappa}_{r}^{\prime}}\text{/}2}}};$ A_(r) = 0; B_(r) = 0; N_(r) = 0; Δκ_(r) = p_(dlpf)Δκ_(r) + (1 − p_(dlpf))Δκ^(′) _(r);  } (138) }

Using the delta scaling factors calculated in (137) and (148), the individual scaling factors of the four wheels could be assigned as in the following if only the smaller wheel is considered

if (Δk_(f) ≧ 0 and Δk_(r) ≧ 0 ) { k₁₀ = 1 − Δk_(f);k₁₁ = 1; k₁₂ = 1 − Δk_(r);k₁₃ = 1; ireff = 1;irefr = 3; } else if (Δk_(f) ≦ 0 and Δk_(r) ≦ 0 ) { k₁₀ = 1 + Δk_(f);k₁₁ = 1; k₁₂ = 1;k₁₃ = 1 + Δk_(r); ireff = 0;irefr = 2; } else if (Δk_(f) ≦ 0 and Δk_(r) ≧ 0 ) { k₁₀ = 1 + Δk_(f);k₁₁ = 1; k₁₂ = 1 − Δk_(r);k₁₃ = 1; ireff = 0;irefr = 3; } else (139) { k₁₀ = 1 − Δk_(f);k₁₁ = 1; k₁₂ = 1;k₁₃ = 1 + Δk_(r); ireff = 1;irefr = 2; }

The least square algorithm for computing Δκ_(f2r) in (134) can be expressed as in the following

if (front and rear reference wheels screening condition holds) { n_(f2r) = n_(f2r) + 1; if (n_(f2r) < p_(totlsstesps)) { A_(f2r) = A_(f2r) + (v_(ireff) + v_(irefr))²/4; B_(f2r) = B_(f2r) + (v_(ireff) ² − v_(irefr) ²)/2;  } else { ${{\Delta \; \kappa_{f2r}} = {\max \; \left( {{- p_{deviation}},{\min \; \left( {p_{deviation},\frac{B_{f2r}}{A_{f2r}}} \right)}} \right)}};$ ${{\Delta\kappa}_{f2r}^{\prime} = \frac{{\Delta\kappa}_{f2r}^{\prime}}{1 - {{{\Delta\kappa}_{f2r}^{\prime}}\text{/}2}}};$ A_(f2r) = 0; B_(f2r) = 0; N_(f2r) = 0; Δκ_(f2r) = p_(dlpf)Δκ_(f2r) + (1 − p_(dlpf))Δκ^(′) _(f2r);  } (140) }

Using the delta scaling factor calculated in (140) the individual scaling factors of the four wheels due to the mismatch between the front reference wheel and the rear reference wheel could be assigned as in the following if only the smaller wheel is considered

if (Δk_(f2r) >= 0) { k₂₀ = 1−Δk_(f2r);k₂₁ = 1−Δk_(f2r); k₂₂ = 1;k₂₃ = 1; iref = irefr; } (141) else { k₂₀ = 1;k₂₁ = 1−Δk_(f2r); k₂₂ = 1 + Δk_(f2r);k₂₃ = 1 + Δk_(f2r); iref = ireff; }

After the adjustment, all the four wheels can have the exact same wheel speeds in a straight driving, however, they could be different from the true wheel rolling radius since the balances are conducted with respect to the reference wheels. The reference wheel is denoted as iref, which could be either the front reference wheel ireff or the rear reference wheel irefr.

If the sensor based longitudinal acceleration is {dot over (v)}_(x-MRF) the following formula

{dot over (v)}_(x-MRF) =a _(x-MRF)+ω_(z-MRF) v _(y-MRF)ω_(y-MRF) v _(z-MRF) +g sin θ_(y-MRF)   (142)

and the time derivative of v_(iref) as {dot over (v)}_(iref), then (126) can be expressed as in the following

$\begin{matrix} {\kappa_{g} = \frac{a_{x - {MRF}} + {\omega_{z - {MRF}}v_{y - {MRF}}} - {\omega_{y - {MRF}}v_{z - {MRF}}} + {g\; \sin \; \theta_{y - {MRF}}}}{{\overset{.}{v}}_{iref}}} & (143) \end{matrix}$

A conditional least square parameter identification method for (143) can be expressed as in the following

if (screening condition holds) { n_(g) = n_(g) + 1;  if (n_(g) < p_(totlsstesps))  { A_(g) = A_(g) + v² _(x−MRF); B_(g) = B_(f) + v_(x−MRF)v_(iref); } (144)  else  { k_(g) = max(−p_(deviation),min(p_(deviation),B_(g) / A_(g))); A_(g) = 0;B_(g) = 0;N_(g) = 0; Δk_(g) = p_(dlpf)Δk_(f) + (1− p_(dlpf))Δk_(g); } }

Itf Δκ_(g) is available, then the final static scaling factors can be calculated as in the following

κ_(stati)=κ_(1i)κ_(2i)κ_(g), i=0,1,2,3   (145)

During aggressive maneuver, for example in an RSC event, the vehicle loading is heavily biased towards the outside wheels, especially during one or two wheel lifting cases. In those cases, there are significant normal loading increases on outside wheels (up to 100% increase in double wheel lift case). Let the tire vertical stiffness be K_(f), the vehicle total mass be M_(t), the distance of the c.g. of the vehicle with respect to the front axle be b_(f) and with respect to front axle be b_(r), wheel base be b, then the vertical compression of the front and rear wheels for double wheel lift case can be calculated as

$\begin{matrix} {{{\Delta \; r_{f}} = \frac{b_{r}M_{t}g}{{bK}_{1}}}{{\Delta \; r_{r}} = \frac{b_{r}M_{t}g}{{bK}_{1}}}} & (146) \end{matrix}$

and in normal driving conditions, they assume the values

$\begin{matrix} {{{\Delta \; r_{f\; 0}} = \frac{b_{r}M_{t}g}{2{bK}_{t}}}{{\Delta \; r_{r\; 0}} = \frac{b_{r}M_{t}g}{2{bK}_{t}}}} & (147) \end{matrix}$

Hence the incremental vertical tire compression changes are the difference between (146) and (147). If taking some common values such as

$\begin{matrix} {{M_{t} = {2704\mspace{14mu} {kg}}},{K_{t} = {300000\frac{N}{m}}},{b_{f} = {1.451\mspace{14mu} m}},{b_{r} = {1.573\mspace{14mu} m}},{b = {3.024\mspace{14mu} m}},{then}} & (148) \\ {{{{\Delta \; r_{f}} - {\Delta \; r_{f\; 0}}} = {\frac{b_{r}M_{t}g}{2{bK}_{t}} = {{0.021\mspace{14mu} m} = {5.57\% \mspace{14mu} {of}\mspace{14mu} R_{f}}}}}{{{\Delta \; r_{r}} - {\Delta \; r_{r\; 0}}} = {\frac{b_{r}M_{t}g}{2{bK}_{t}} = {{0.023\mspace{14mu} m} = {6.10\% \mspace{14mu} {of}\mspace{14mu} R_{r}}}}}} & (149) \end{matrix}$

where R_(f)=R_(r)=0.377 m denote the nominal rolling radii of the tires. Hence for double wheel lift case, the outside tires could have around 5% reduction in rolling radius. For this reason, a scaling factor is computed to compensate this

if (|θ_(xr) |≦ p_(35th)p_(rollgrad)) { κ_(dyn0) = 1; κ_(dyn1) = 1; κ_(dyn2) = 1; κ_(dyn3) = 1;   } else { if (θ_(xr) < 0) (150) { ${\kappa_{{dyn}\; 0} = {\max \; \left( {{1 - p_{deltalatamax}},{1 - {\frac{{\theta_{xr}} - {p_{35{th}}p_{rollgrad}}}{p_{40{th}}p_{rollgrad}}p_{deltalatamax}}}} \right)}};$ κ_(dyn2) = κ_(dyn0); κ_(dyn3) = 1; κ_(dyn3) = 1; else { ${\kappa_{{dyn}\; 1} = {\max \; \left( {{1 - p_{deltalatamax}},{1 - {\frac{{\theta_{xr}} - {p_{35{th}}p_{rollgrad}}}{p_{40{th}}p_{rollgrad}}p_{deltalatamax}}}} \right)}};$ κ_(dyn3) = κ_(dyn1); κ_(dyn0) = 1; κ_(dyn2) = 1; } where κ_(dyni) (for i=0,1,2,3) denotes the dynamic scaling factor for the ith wheel.

The final individual scaling factors can be calculated as in the following

κ_(i)=κ_(stati)κ_(dyni)=κ_(1i)κ_(2i)κ_(g)κ_(dyni), i=0,1,2,3   (151)

In the following example the above computations can be used for compensating the wheel speed sensor signals. The front-left wheel and the rear left wheels are wheels presumed to have smaller rolling radii, then the compensated wheel speeds using the above two stage compensations can be expressed as in the following for Δκ_(f2r)<0 case

w ₀ ^(comp)=(1−|Δκ_(f)|)(1−|Δκ_(f2r)|)w ₀

w ₁ ^(comp)=(1−|Δκ_(f2r)|)w _(i)

w ₂ ^(comp)=(1−|Δκ_(r)|)w ₂

w₃ ^(comp)=w₃   (152)

and for Δκ_(f2r)≧0 case

w ₀ ^(comp)=(1−|Δκ_(f)|)w ₀

w₁ ^(comp)−w₁

w ₂ ^(comp)=(1−|Δκ_(r)|)(1−|Δκ_(f2r)|)w ₂

w ₃ ^(comp)=(1−Δκ_(f2r)|)w ₃   (153)

It can be seen that in either case, there is one-wheel which does not need to be compensated (rear-right wheel for (151) and front right wheel for (152)). If considering κ_(g) and κ_(dyni), further compensation for (160) may be achieved as in the following

w ₀ ^(comp)κ_(g)(1−|Δκ_(f)|)(1−|Δκ_(f2r)|)κ_(dyn0) w ₀=κ₀

w ₁ ^(comp)=κ_(g)(1−|Δκ_(f2r)|)κ_(dyn1) w ₁=κ₁ w ₁

w ₂ ^(comp)=κ_(g)(1−|Δκ_(r)|)κ_(dyn2) w ₂=κ₂ w ₂

w₃ ^(comp)=κ_(g)κ_(dyn3)w₃=κ₃w₃   (154)

or for (153) as in the following

w ₀ ^(comp)=κ_(g)(−|Δκ_(f)|)κ_(dyn0) w ₀=κ₀ w ₀

w₁ ^(comp)=κ_(g)κ_(dyn1)w₁κ₁w₁

w ₂ ^(comp)=κ_(g)(−|Δκ_(r)|)(1−|Δκ_(f2r)|)κ_(dyn2) w ₂=κ₂

w ₃ ^(comp)=κ_(g)(1−|Δκ_(f2r)|)κ_(dyn3) w ₃=κ₃ w ₃   (155)

The detailed algorithm together with the screening conditions used in the algorithms may be found in the detailed logic section.

Referring now also to FIGS. 12A, 12B, and 12C, one advantage of the present invention is illustrated. As mentioned above, the longitudinal velocity, the lateral velocity determined above may be used to determine the sideslip of the vehicle. In FIG. 12A, a base vehicle without stability control is illustrated. As can be seen, before achieving stability the vehicle may completely slide sideways. In a traditional stability control system shown in 12B, the vehicle may slide sideways and leave the roadway before regaining control. This is because traditional stability control may be too late and the vehicle may overshoot the desired course. Using the present invention illustrated in 12C, an advance yaw stability control using precise sideslip information reacts quicker than standard stability control so that the desired course or lane may be maintained.

The vehicle Sideslip Angle at the Rear Axle

As mentioned above, the sideslip angle β_(LSC) of the vehicle body and its time derivative are used in LSC. β_(LSC) is the sideslip angle β_(mrfra) of the vehicle body at its real axle location, measured along the moving road frame, measures the lateral stability of a moving vehicle. The aforementioned sideslip angle can be defined as in the following

$\begin{matrix} {\beta_{LSC} = {\tan \left( \frac{v_{mrfray}}{v_{mrfrax}} \right)}} & (156) \end{matrix}$

where v_(mrfraz) and v_(mrfrax) are the longitudinal and lateral velocities of the vehicle body at its rear axle measured along the x-axis and the y-axis of the moving road frame.

β_(LSC) thus defined can be related to the signals calculated from the sensor measurements such as the lateral acceleration, the yaw rate, the bank and slope angles of the moving road frame and the reference velocity as in the following:

$\begin{matrix} {{d\; \beta_{{LSC}\; 2}} = {\left\lbrack {{- \omega_{mrfz}} + \frac{{a_{mrfyra}v_{mrfrax}} - {g\; \sin \; \theta_{mrfx}\cos \; \theta_{mrfy}}}{v_{mrfrax}}} \right\rbrack {\quad{{\left\lbrack {1 + {\cos \left( \beta_{{LSC}\; 2} \right)}} \right\rbrack - {\frac{{\overset{.}{v}}_{mrfrax}}{v_{mrfrax}}{\sin \left( \beta_{{LSC}\; 2} \right)}\beta_{LSC}}} = \frac{\beta_{{LSC}\; 2}}{2}}}}} & (157) \end{matrix}$

Both (156) and (157) can be used to compute β_(LSC). (146) requires the computation of the lateral velocity while (157) does not need to have a separate computation of lateral velocity. In the following equation (157) is used to compute β_(LSC). Due to the inevitable sensor uncertainties in lateral acceleration and yaw rate sensors, and errors in the computed road bank and slope, a pure integration of (157) would generate large error.

In step 924 a dynamic sideslip angle may be determined, such a dynamic computation calculate the integration of (157) with passing through a high-pass-filter. Such integration plus high pass filter is called an anti-integration-drift (AID) filter. Such computation can be done as in the following

$\begin{matrix} {{\left. {{\beta_{{LSC}\; 2}\left( {k + 1} \right)} = {{d_{1}{\beta_{{LSC}\; 2}(k)}} - {d_{2}{\beta_{{LSC}\; 2}\left( {k - 1} \right)}}}} \right\rbrack + {n_{1}\left\lbrack {{d\; {\beta_{{LSC}\; 2}\left( {k + 1} \right)}} - {d\; {\beta_{{LSC}\; 2}\left( {k - 1} \right)}}} \right\rbrack}}{{\beta_{{LSC}\; 2}(0)} = {2{\beta_{mrflinr}(0)}}}{{\beta_{{LSC} - {DYN}}\left( {k + 1} \right)} = \frac{\beta_{{LSC}\; 2}\left( {k + 1} \right)}{2}}} & (158) \end{matrix}$

where d₁,d₂ and n₁ are the coefficients of the AID filter.

The computation in (158) captures the dynamic portion of the sideslip angle, which could be accurate during aggressive maneuvers but it could be either under-estimated or over-estimated in less dynamic but unstable situations such as driving on a low mu road surface. For this reason, (158) above does not have much use unless a compensation for less dynamic portion can be calculated. In order to conduct less dynamic compensation, the similar scheme used in the RSG unit of the ISS may be conducted.

The less-dynamic portion might calculated from the previous computation of v_(mrfray) and v_(mrfraz) through the following steady-state-recovery (SSR) filtering

$\begin{matrix} {\left. {{\beta_{{LSC} - {SS}}\left( {k + 1} \right)} = {{d_{1}{\beta_{{LSC} - {SS}}(k)}} - {d_{2}{\beta_{{LSC} - {SS}}\left( {k - 1} \right)}}}} \right\rbrack + {{f_{1}\left\lbrack \frac{v_{mrfray}}{v_{mrfrax}} \right\rbrack}\left( {k + 1} \right)} + {{f_{2}\left\lbrack \frac{v_{mrfray}}{v_{mrfrax}} \right\rbrack}(k)} + {{f_{3}\left\lbrack \frac{v_{mrfray}}{v_{mrfrax}} \right\rbrack}\left( {k - 1} \right)}} & (159) \end{matrix}$

with proper filter coefficients f₁,f₂ and f₃.

In step 928 the final sideslip angle may be expressed as

β_(LSC)=β_(LSC-DYN)+β_(LSC-SS)   (160)

Advantageously, the present invention provides an improved sideslip angle determination by fully utilizing the available measured and calculated signals from the ISS system. Such computation corrects for gravity and nonlinearities as well as steady state sideslip angle factors and removing low frequency drift due to the inevitable integration error due to sensor signal uncertainties such as zero-offset and temperature dependent offset drifting, error in the computed signals. Most importantly the sideslip angle calculated here is defined in the moving road plane, which reflect the true objective of the LSC and which is different from many existing calculations that are defined along the body-fixed frame.

In step 930 the vehicle and the various control systems described above may then be controlled in response to the final sideslip angle. As mentioned above, engine, brakes, steering and other actuators may be controlled to change the dynamics of the vehicle.

While particular embodiments of the invention have been shown and described, numerous variations and alternate embodiments will occur to those skilled in the art. Accordingly, it is intended that the invention be limited only in terms of the appended claims. 

1. A method of controlling a vehicle comprising: determining a front lateral tire force; determining a rear lateral tire force; determining a linear sideslip angle from the front lateral tire force and the rear lateral tire force; determining a linear lateral velocity in response to the linear sideslip angle; and controlling the vehicle in response to the linear lateral velocity.
 2. A method as recited in claim 1 further comprising determining a longitudinal kinematics constraint and determining the linear lateral velocity constraint in response to the linear sideslip angle and the longitudinal kinematics constraint.
 3. A method as recited in claim 2 further comprising determining the longitudinal kinematics constraint in response to a longitudinal velocity, a yaw rate and a longitudinal acceleration.
 4. A method as recited in claim 1 further comprising determining a road constraint alignment and determining the linear lateral velocity constraint in response to the linear sideslip angle and the road constraint alignment.
 5. A method as recited in claim 4 further comprising determining the road constraint alignment in response to a longitudinal velocity, a pitch rate and a roll rate.
 6. A method as recited in claim 1 further comprising determining a lateral velocity error due to longitudinal velocity error.
 7. A method as recited in claim 1 further comprising determining a lateral velocity error due to an augmented acceleration.
 8. A method as recited in claim 7 wherein the augmented acceleration comprises a lateral augmented acceleration error.
 9. A method as recited in claim 8 wherein the lateral augmented acceleration error is determined in response to an attitude error.
 10. A method as recited in claim 9 wherein the attitude error comprises a pitch angle error.
 11. A method as recited in claim 9 wherein the attitude error comprises a roll angle error.
 12. A method as recited in claim 7 wherein the augmented acceleration comprises a longitudinal augmented acceleration error.
 13. A method as recited in claim 12 wherein the longitudinal augmented acceleration error is determined in response to an attitude error.
 14. A method as recited in claim 13 wherein the attitude error comprises a pitch angle error.
 15. A method as recited in claim 7 wherein the augmented acceleration comprises a vertical augmented acceleration error.
 16. A method as recited in claim 15 wherein the vertical augmented acceleration error is determined in response to an attitude error.
 17. A method as recited in claim 16 wherein the attitude error comprises a pitch angle error.
 18. A method as recited in claim 16 wherein the attitude error comprises a roll angle error
 19. A method of controlling a vehicle comprising: determining a yaw rate from a yaw rate sensor; determining a pitch rate from pitch rate sensor; determining a roll rate from a roll rate sensor; determining a pitch rate from a pitch rate sensor; determining a lateral acceleration from a lateral acceleration sensor; a vertical acceleration from a vertical acceleration sensor; determining a front lateral tire force in response to the yaw rate, pitch rate, roll rate, lateral acceleration and vertical acceleration; determining a rear lateral tire force in response to the yaw rate, pitch rate, roll rate, lateral acceleration and vertical acceleration; determining a linear sideslip angle from the front lateral tire force and the rear lateral tire force; determining a linear lateral velocity in response to the linear sideslip angle; and controlling the vehicle in response to the linear lateral velocity.
 20. A method as recited in claim 19 further comprising determining a roll inertia and wherein determining a front lateral tire force comprises determining the front lateral tire force in response to the yaw rate, pitch rate, roll rate, lateral acceleration, vertical acceleration and roll inertia.
 21. A method as recited in claim 19 further comprising determining a roll inertia and wherein determining a rear lateral tire force comprises determining the rear lateral tire force in response to the yaw rate, pitch rate, roll rate, lateral acceleration, vertical acceleration and roll inertia.
 22. A method as recited in claim 19 further comprising determining a pitch inertia and wherein determining a front lateral tire force comprises determining the front lateral tire force in response to the yaw rate, pitch rate, roll rate, lateral acceleration, vertical acceleration and pitch inertia.
 23. A method as recited in claim 19 further comprising determining a pitch inertia and wherein determining a rear lateral tire force comprises determining the rear lateral tire force in response to the yaw rate, pitch rate, roll rate, lateral acceleration, vertical acceleration and pitch inertia 